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ABSTRACT 


The design, development, fabrication and testing of a siringle-type terrestrial solar 

rt 2 • • 

cell module which produces 98 watts /m of exposed module area at 1 kW/m insolation 
and 61®C are reported. These modules make it possible to easily incorporate- photovoltaic 
power generation into the sloping roofs of residential or commercial buildings by simply 
nailing the modules to the. plywood roof sheathing. 

This design consists of nineteen series-connected 53 mm diameter solar cells arranged 
in a closely packed hexagon configuration. These cells are individually bonded to the 
embossed surface of a 3 mm thick thermally tempered hexagon-shaped piece of ASG 
SUNADEX glass. Monsanto SAFLEX polyvinyl butyral is used as the laminating adhesive. 
RTVII functions as the encapsulant between the underside of the glass superstrate 
and a rear protective sheet of 0.8 mm thick TEXTOLITE. The semi— flexible portion of 
each shingle module is a composite laminate construction consisting of outer layers of 
B. F. Goodrich FLEX SEAL and an epichlorohydrin closed cell foam core. 

The module design has satisfactorily survived the JPL - defined qualification testing 
program which includes 50 thermal cycles between -40 and +90°C, a seven day tempera- 
ture - humidity exposure test and a mechanical integrity test consisting of a bidirectional 
cyclic loading at 2390 Pa (50 Ib/ft^) which is intended to simulate loads due to a 45 m/s 
flOO mnhl wind. 
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SECTION 1 
SUMMARY 


The shingle-type solar cell module shown in Figure 1-1 has been designed and developed 
during the one year period covered by this final report. A total of 50 modules have 
been fabricated for delivery to JPL. Six of these have been subjected to an environ- 
mental testing program consisting of 50 thermal cycles between -40 to +90°C , a one 

week temperature-humidity exposure, and a 100 cycle bidirectional mechanical integrity 

2 

test to simulate snow and wind loading at 2390 Pa (50 lb /ft ). This module design sand- 
wiches the interconnected solar cells between a sheet of tempered SUNADEX glass on 
the front surface and a sheet of fiberglass /epoxy on the rear side. The 19 series- 
connected solar cells are bonded to the embossed surface of this glass with polyvinyl 
butyral (PVB) film and the space between the covers is filled with RTVII which functions 
as the primary encapsulant. The semi- flexible portion of each module is a laminate 
constructed of B. F. Goodrich FLEXSEAL outer skins and an epichlorohydrin closed 
cell foam core. A two-sided printed wiring board, which is sandwiched within this 
laminate, provides the connection between the ends of the solar cell circuit string and 
the four output terminals of the module. These terminals overlap, negative on positive, 
and are interconnected with a machine screw /flat washer to produce the shingle installa- 
tion shown in Figure 1-2. 

9 

An average module electrical performance of 5.79 watts was measured at 100 mW/cm 
insolation and at 28°C using JPL-supplied Terrestrial Secondary Standard No. 025 as 
the reference. This represents an overall module efficiency of 11.4 percent based on 
the exposed module area. If this performance is reflected to the calculated Nominal 
Operating Cell Temperature (NOCT) of 61°C, the resulting module output yields an 
areal specific power output of 98 watts/m of exposed module area at 100 mW/cm'^ 
insolation and at the NOCT. 
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(a) Front (Photo No. VF 78-278A) 



(b) Back (VF Photo No. 77-557A) 

Figure 1-1. Photographs of Shingle Solar Cell Module 
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Figure 1-2. Arrangement of Modules on Simulated Roof Structure 

(Photo No. VF 78-237B) 



SECTION 2 


INTRODUCTION 



SECTION 2 
INTRODUCTION 

The general scope of work under this contract involves the design , development , fabri- 
cation and testing of a solar cell module which is suitable for use in place of shingles 
on the sloping roofs of residential or commercial braidings. Modules of this type em- 
ploy a semi-flexible substrate which is suitable for mounting on an independent rigid 

surface such as plywood roof sheathing. As specified in the contract statement of work, 

2 

these modules shall be capable of producing an electrical power output of 80 W/m of 

o 2 

installed module area at a module temperature of 60°C with an insolation of 1 kW/m . 

The installed weight of these shingle-type modules shall not exceed 250 kg/kW of peak 
power output. As a design goal these modules shall be designed for a service life of 
at least 15 years. An implicit requirement is that the shingle not sustain damage during 
the normal handling associated with installation on a roof. The vulnerability to the 
localized bearing loads associated with walking or kneeling on the installed shingles 
does not constitute a design requirement but will be assessed as part of this develop- 
ment effort . The program is organized into seven major tasks as given below . 

Description 

Substrate Evaluation and Testing 
Solar Cell Tray Evaluation and Testing 
Module Interconnection and Testing 
Shingle Module Design 

Fabrication and Acceptance Testing of Modules 

Qualification Testing of Modules 

Analytical Model of a Zero Depth Solar Photovoltaic 
Concentrator 

During the initial phases of the program, the activities on Tasks 1, 2 and 3 involved 
the investigation of a variety of shingle module implementation approaches. 

Considerable effort was expended in an attempt to develop an approach which used 
methyl methacrylate (MMA) as the sole solar cell encapsulant. This proved to be 


Task-No. 

1 

2 

3 

4 

5 

6 
7 
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impractical because the relatively high modulus MMA could not elastically accommodate 
the strains at the specified low temperature extreme of -40°C. Attempts to buffer 
the solar cells with a transparent silicone conformal coating prior to embedding within 
MMA also proved to be unsuccessful. 

This activity, as well as the evaluation of various module-to-module interconnection 
concepts , has been reported in the first two quarterly reports published under this 
contract. These results will not be repeated here, but rather the emphasis of this 
final report will be on the design, fabrication, and testing of the tempered glass covered 
shingle module. 

Task 7 was added after the completion of the initial contract effort. The results of 
this additional task activity, which involved the development of an analjrtical model of 
a zero depth solar photovoltaic concentrator, are reported in Appendix A to this docu- 
ment. 

Fifty shingle solar cell modules were deKvered to JPL on May 19, 1978. Forty of these 
modules were mounted on a simulated roof structure and have been undergoing out- 
door exposure testing at the JPL Field Test Site since October 17, 1978. 
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SECTION 3 

TECHNICAL DISCUSSION 

3.1 DESCRIPTION OF SELECTED DESIGN 

3.1.1 GENERAL DESCRIPTION 

The selected design for the shingle solar cell module is represented by the assembly 
drawing shown in Figure 3-1. This module consists of two basic functional parts; an 
exposed rigid portion which contains the solar cell assembly, and a flexible portion which 
is overlapped by the higher courses of the roof installation. The design of the shingle 
module provides a closely-packed array of 19 series-connected solar cells. A minimum 
separation of 0.5 mm (0.020 inch) is maintained between adjacent ceRs by assembly tool- 
ing which positions the cells prior to bonding to the glass coverplate. The position of the 
four output terminals of the module has been established to permit the connection of the 
negative terminals of one course on the roof directly to the positive terminals of the 
course below. The method of connection, which uses a. machine screw and flat washer, 
is discussed in Section 3.1.6. 

As shown in Sections C-C .and D-D of Figure 3-1, the top substrate FLEXSEAL skin over- 
laps , and is bonded to , -the glass coverplate to form a weather-tight joint around the 
upper three sides of the hexagon. The bottom skin and printed wiring board ave sand- 
wiched between the glass coverplate and the bottom fiberglass /epoxy cover to produce, 
a similar seal around these three edges on the bottom. The exposed edges of the glass 
coverplate are sealed with a bead of RTV102 which is applied between the coverplate and 
the bottom fiberglass /epoxy cover as. shown in Section F-F. 

The key features of this shingle module design are summarized in Table ‘ 3-1. The cal- 
culated module output of 4.95 watts at the Standard Operating- Conditions (SOC) which 
' include a calculated NOCT of 61°C, yields an areal specific output of 98 watts/m^ of ex- 
posed module area. This is 23 percent better than the minimum requirement of the con- 
tract. The specific weight of the module is 202 kg/kW of peak power output at SOC as 
compared to a maximum specified value of 250 kg/kW. 
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Table 3-1. Key Featiires of Shingle Module Design 


Parameter 

Value 

Total Solar Cell Area 

419.2 cm^ 

Exposed Module Area 

507.0 cm^ 

Packing Factor 

0.827 

Electrical Power Output at the Maximum 
Power Point 


At 1 kW/m^ and 28°C 

5.79 Watts* 

. At Standard Operating Conditions 

4.95 Watts 

Module Weight 

1.00 kg* 


■ * Average of 50 modules delivered 


3.1.2 SUBSTRATE CONFIGURATION 

The semi-flexible substrate portion of the shingle is of laminar construction as shown 
in Figure 3-2. The two outer-sMns of this substrate are FLEXSEAL polyester scrim re- 
inforced HYPALON . This material is white in color and provides the weather-resistance 
properties required to meet the 15 year service life goals. A center core of closed cell 
epichldrohydrin foam (Rubatex No. R-^473-E) provides a low-density, high-temperature 
resistant filler material to achieve a nearly uniform thickness of the entire surface area 
of the shingle. The substrate also affords protection to the flexible printed wiring board 
which is sandwiched between the bottom sMn and the core. This double sided printed 
wiring board, which carries both the positive and negative terminations for the module, 
is made from material identified by GE designation FLGF 0.006 C 2/2 which consists of 
2 oz’/ft2 copper foil on both sides of a 0,006. inch thick fiberglass /epoxy substrate. This 
copper is etched away to form the bus strip patterns. The calculated series resistance 
of this copper bus network is the 4.6 m at 20°C, which represents a negligible- power 
loss. 

The entire composite substrate is bonded together with the B.F. Goodrich adhesive sys- 
tem described in Table 3-2. This contact adhesive and associated primers are tradition- 
ally used as part of the FLEXSEAL roofing system. Two different primers have been 
specified depending upon the nature of the surface to be bonded. In either case the 
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primer is applied in a very thin coating and allowed to dry thoroughly prior to the appli- 
cation of the contact adhesive to both surfaces to be joined. The edge sealer (A1436-B) 
is applied to the outer edges of the substrate to seal the exposed foam core material. 



Figure 3-2. Section Through Shingle Module Substrate 


Table 3-2. Substrate Adhesive System 


Identification No. * 

Description /AppUeation 

A 1104-B 

Primer for Non-porous Surfaces 

A 178-B 

Primer for Porous Surfaces 

CA-1056 

Contact Adhesive (all surfaces) 

A 1436-B 

Edge Sealer 


* B. F. Goodrich Co. 

General Products Division 
Solon, Ohio 44139 
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A load-deflection test was performed on, a segment of the substrate consisting of a top 
and bottom sMn of polyester scrim rranforced FLEXSEAL- bonded to a cord of epicMoro* 
hydrin foam. The resulting load-deflection curve with a 14.3 mm (0.563 inch) diameter 
bearing surface is as given in Figure 3-3, The use of these data in the determination of 
module-to-module interconnection joint contact force will be discussed later. 

3.1.3 SOLAR CELL SELECTION 

A Spectrolab solar cell was. selected for- use in this module design on the basic of the 
lowest sp.ecifie cost of the delivered power output. The distribution of the electrical out- 
put of the cells -received from Spectrolab is given in Table 3-3. As measured by Spectro- 
lab these cells have an average electrical output of 569 mA at 0.‘475 volts when measured 
at 1 kVf/m^ insolation and 38®C. A random sample of 66 ceUs from 'this group yielded an 
average performance of 595 mA at 0.475 volts under these same test conditions when mea- 
sured by GE using the Large Area Pulse Solar Simulator (LAPSS) with JPL - supplied 
Terrestrial Secondary Standard No. 025 as the reference. 

Scanning electron micrographs of the front and rear contacts of these cells are shown 
in Figures 3-4 and 3-5, respectively. These micrographs clearly reveal the nature of 
the printed contact surface. 

3.1.4 SOLAR CELL INTERCONNECTOR 

The solar cell interconnector shown in Figure 3-6 is fabricated from nominal SOfim (0.002 
inch) thick soft copper foil (AUoy No. 110) which is subsequently solder plated to a 
thickness of 13 to 18 ^m (0.0005 to 0,0007 inch) on both surfaces. The resistance of 
this interconnector measured between the two ”N” joints and the three "P” joints is 1.67 
mQ at 25°C. At 60°C-this series resistance loss amounts to 0.3. percent of the cell maxi- 
mum power output. Figure 3-7 shows a typical ”N” contact solder- joint which was made 
by reflowing the solder platingcon the interconnector with a hand soldering iron. 

3.1.5 MODULE ENCAPUSLATION 

The details of the encapsulation surrounding the solar cell assembly are shown in Figure 
3-8. The solar cells are individually bonded to. the. underside of the glass. coverplate 
with disks of Monsanto SAFLEX SR-10 PVB fidm. This glass coverplate is fabricated from 
3 mm thick ASG SUNADEX low-iron glass (0.01% iron-roxide content) which is cut to the 
hexagon shape and then thermally tempered to provide the flexual strength required to 
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Figure 3-3. Local Bearing Load-Deflection Curve for Substrate 



Table 3-3. Electrical Performance of Production Solar Cells 
as Measured by Spectrolab 


Range of 
Output Current 
Group (mA at 0.475 V) 

Quantity of 
Solar CeUs 
In Group 

510 - 519 

69 

520 - 529 

95 

530 - 539 

108 

540 - 549 

115 

550 - 559 

.120 

560 - 569 

.127 

570 - 579 

164 

580 - 589 

147 

590 - 599 

146 

600 - 609 

97 

610 - 619 

72 

620 - 629 

27 

630 - 639 
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sustain the bearing loads associated with walking or kneeling. The transmission of 
this glass is compared with that of ’ASG LO-IRON (0.05% iron-oxide content) in Figure 
3-9. These data do not indicate a clear transmittance advantage associated with the use 
of the SUNADEX glass, but this selection was made because of the embossed surface 
texture of this glass and its influence on the enhanced output of the module as discussed' 
in Section 3.4. 


The primary encapsulation around the solar cells is provided by RT'V 11, which is a 
white pourable conformal coating. This dimethyl silicone compound fills the space sur- 
rounding the solar cells and interconnectors and bonds the front glass coverplate to the 
rear protective sheet of fiberglass /epoxy. This rear sheet is required to prevent damage 
to the module from sharp objects such as nails which are everpresent during the instal- 
lation of a shingle roof. The adhesion of the RT'V 11 to both the glass coverplate and the 
rear TEXTOLITE sh6et is increased by the use of GE ,primer number SS-4044 followed by 
a thin coating of RTV 108 (10 to 20’ percent by weight mixed with heptane). 

3.1.6 MODULE-TO-MODULE INTERCONNECTION 

The mod\ile-to-module interconnection is accomplished as shown in Figure 3-10. This 
basic concept relies on the development of high contact pre^ssure under three conical pro- 
jections , which .are part of solder plated copper bosses within each mating shingle , to 
achieve a low-resistance, environmental stable connection. 
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(a) Front Contact Interface (b) Front Contact 

With Active Area 

Figure 3-4. Scanning Electron Micrographs of Solar Cell Front Contact 



Figure 3-5. Scanning Electron Micrograph of Solar Cell Back Contact 
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Figure 3-6. Solar Cell Interconnector 



Figure 3-7. Typical "N” Contact Solder Joint 



Figure 3-8. Module Encapsxilation 





Figure 3-9. Transmission of ASG Low-Iron Soda-Lime Glasses 
Without Anti-Reflection Coating 
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PRINTED 
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Figure 3-10. Module- to- Module Interconnection 
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The two positive terminals of each shingle are provided with these copper bosses which 
are soldered to the circular pads on the top of the printed circuit board. A 6-32 
TENNUT, which has been modified by shortening the barrel and removing the prongs, 
is centered within each of these positive terminals with an insulating sleeve. Both this 
sleeve and the TEENUT are made captive by bonding to each other and to the copper 
boss. The TEENUT is thus insulated from the positive terminal by the sleeve and by 
the bottom FLEXSEAL skin under the circvut board. The two negative terminals of 
each module are provided with thinner copper bosses which are soldered to the bottom 
of the printed circuit board at the circular pad mounting locations. Upon assembly, con- 
tact pressure is developed between these two overlapping copper bosses by compressing 
the substrate core with a 6-32 flat head machine screw and 16.5 mm (0.650 inch) dia- 
meter flat washer. The screw is driven into the TEENUT until the top surface of the 
washer is flush with the top surface of the outer FLEXSEAL skin as shown in Figure 
3-11. Under these conditions the 1.8 mm (0.070 inch) thick washer will have developed 
a total compressive force of 52 N (see Figure 3-4) or 17.3 N per conical projection. 

This is adequate force to provide a reliable high-pressure, gas-tight contact. The in- 
sulating sleeve around the screw thread prevents electrical contact with the negative 



Figure 3-11. Installation of Module-to-Module Interconnector 
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terminal of the module so that the exposed screw head /washer /TEEN UT are electricaUy 
neutral after installation. This interconnector design concept is also capable of accom- 
modating up to 1.8 mm (0.07 inch) of misregistration between the centerlines of over- 
lapping terminal bosses. 

3.2 SYSTEM INSTALLATION CONSIDERATIONS 

The shingle modules described above are installed by overlapping in the manner shown 
in Figures 3-12 and 3-13. The four electrical terminals of each module are interconnected 
as previously described using a flat-head machine screw and washer to apply the contact 
force. The shingles are attached to the roof sheathing by nailing through the substrate 
at two places per shingle with ordinary roofing nails as shown in Figure 3-14. The 
module- to-module interconnectors between overlapping layers form a series /parallel 
matrix of interconnected modules, as shown in Figure 3-15, with the current increasing 
as modules are added in the parallel direction across the length of the roof from gable- 
to- gable and voltage increasing as modules are added in the series direction along the 
slant height of the roof from eave-to- ridge. At both the positive and negative termin- 
ations of each solar cell circuit it is necessary to attach conductors to carry the current 



Figure 3-12. Typical Overlap Between Courses 
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MODULE-TO-MODULE 




Figure 3-14. Nailing of Shingle Module to Roof Sheathing 
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Figure 3-15. Electrical Schematic of Module-to-Module Interconnection 


to the central collection point which is at the ridge of the roof. A negative circuit ter- 
mination which is typical of that required at the eave of the roof is shown in Figure 
3-16. The connection between the negative terminals of adjacent modules in the first 
active course of the roof is made by underlaying a special shingle which contains a built- 
in printed circuit board to make this jumper connection. At every fifth parallel-connected 
module a tap is made from this common negative bus by soldering a flat-conductor cable 
to this special shingle and rimning this cable under the shingles to the ridge of the roof 
as shown in Figure 3-17. In this way it is possible to limit the current density in the 
common negative bus which is the negative conductor within the shingle modules , and 
control the series resistance losses for the overall installation. Similar terminations are 
required at the ridge of the roof and at electrically common bus connections which may 
be required at intermediate points along the slant height because of requirements to 
provide a specified voltage within the physical constraints imposed by roof dimensions. 

In this latter case modules with opposite polarity for the output terminals are used to 
built-up voltage in the opposite sense as the roof is advanced toward the ridge. 
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Figure 3-16. Second Course Containing Negative Circuit Termination 



Figure 3-17. Installation of Flat Conductor Cable 


Of P"®’* 
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The determination of the size for the rectangular roof siirface area required for an in- 
stallation of these shingle modules is based on the general arrangement given in Figure 
3-18. Using this nomenclature, W is the reqxiired length of the roof measured from 
gable- to- gable and S is the slant height meeusured from eave- to- ridge. The values of 
these dimensions are given by: 

W = E (3 Np + 1/2) 
and 

S = E (Ng + 1) cos 30*? + T, 
where: 

E = edge dimension of hexagon (mm) 

= 139.7 mm 

Np = number of parallel-connected modules across the length of the roof 

N_ = number of series-connected modules along the slant height at the roof 

T = space at ridge required to bus together top terminals (mm) 

= 63.5 mm 

The cross hatched area on Figure 3-18 represents wasted space which detracts from the 
overall packing efficiency of the Installation. For a typical installation of 1872 modules 
which are arranged as three circuits, each with 26 series connected modviles and 24 
parallel-connected modviles, the dimension W is 10.128 m. The slant height of the roof 
(S) is 9.621 m since there are three circviits in this direction which yields a value of 
Ng = 78 for the installation. Thus, the required roof area of 97.44 m^ for this typical 
installation results in an overall area utilization of 0.805, which is defined as the ratio 
of solar cell area to rectangfular roof area. In more general terms Figure 3-19 gives this 
area utilization as a function of the number of parallel connected modules and the total 
solar cell area. For smaller instaUations , these results indicate the importance of using 
roof sections which are narrow from gable- to- gable and long from eave-to-ridge. 
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Figure 3-18. Arrangement of Shingle Modules on a Rectangular Roof 


3.3 MODULE FABRICATION 
3.3.1 INTRODUCTION 

The module assembly was begun on February 28, 1978 with the bonding of the 19 cells 
for module Serial No. SM-1 and completed by May 5, 1978 with the encapsulation of 
Serial No. SM-59. During this time several changes in design and process control had 
been instituted to improve the performance and producibility of the shingle modules. 
Table 3-4 summarizes, by serial number, the significant differences of this type which 
are latent within each module. The serial numbers which contain an ”A” have been sub- 
jected to an autoclave processing step as described in Section 3.3.3. As a consequence 
these modules , which have an unstable solution of air within the PVB , are subject to the 
appearance of bubbles as a result of slight elevations in temperature. Such a condition 
did occur on a few of these modviles, notably on serial number SM-IOAQ when subjected 
to the thermal cycling test. As reported in Section 3.5.1, a change in maximum power 
output of only 2 percent resiilted from the appearance of these bubbles. The serial 
numbers which contain a "Q” have been subjected to the qualification testing program 
as reported in Section 3.5. ' 
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Figure 3-19. Area Utilization for Shingle Modules Mounted on a 

Rectangular Roof 
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Table 3-4. Svimmary of Module Fabrication Differences 


Module 

Serial 

Number 

SAFLEX 

Type 

Autoclave 

Cycle 

Primer on 
Rear Surface 
Of Glass 

Rework 
To Replace 
Broken Cells 

SM-IQ 

PT- 

-10 

No 

RTV108 

No 

-2Q 



No 

RTV108 

No 

-3 


i 

No 

RTV108 

Yes (1) 

-40 



No 

RTV108 

No 

-5A 



Yes 

SS4044 -I- RTV108 

Yes (3) 

-6Q 



No 

RTV108 

No 

-7Q 



No 

RTV108 


i 

-lOAQ 



Yes 

RTV108 



-12A 



Yes 

RTV108 



-13A 

PT- 

-10 

Yes 

SS4044 -I- RTV108 



-17A 

SR- 

-10 

Yes 

SS4044 RTV108 



-18 


No 

RTV108 



-19 


No 

RTV108 



-20A 



Yes 

SS4044 + RTV108 



-21 



No 

RTV108 


r 

-22Q 





SS4044 + RTV108 



-23 





j 


No 

-24 





t 

i 

Yes (1) 

-26 







Yes (3) 

-27 







No 

-28 







i 

i 

-29 









-30 









-31 







1 

r 

-32 







No 

-34 







No 

-35 







No 

-36 







Yes (1) 

-37 







No 

-38 







No 

-39 







No 

-41 







Yes (1) 

-42 







Yes (1) 

-43 







No 

-44 







i 

i 

-45 









-46 









-47 







1 

f 

-48 







No 

-49 







Yes (1) 

-50 







No 

-51 







i 

i 

-52 









-53 









-54 









-55 









-56 









-57 









-58 

y 



f 


f 



-59 

SR 

-10 

No 

SS4044 RTV108 

No 
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A change in SAFLEX type from PT-10 to SR- 10 was made with serial number SM-17A. 

The reasons for this change are discussed in Section 3.3.3. The use of a two step 
priming system on the rear side of the glass coverplate was instituted as indicated in 
the table to enhance the adhesion of the RTVll pottant. A rework procedure to replace 
cracked solar ceUs which resulted from the PVB bonding operation was performed on 
ten of the modules as indicated in parentheses . With the exception of serial number 
SM-3 aU of these replacement cells were bonded to the glass coverplate with RTV655. 

In the former case the replacement cell was laminated to the glass coverplate with 
SAFLEX SR- 10 using a vacuvim bag to apply pressure while locally heating the cell by 
passing a forward current of approximately 25 amperes through the cell. This rework 
method proved to be unsatisfactory due to the nonuniform temperature distribution over 
the surface of the cell with hot spots at the ”N” contact solder joints. 

3.3.2 SOLAR CELL SOLDERING 

Soldering of the interconnector to the solar cell ”N" contacts is the first step in the 
assembly sequence. The solder-plated copper interconnector is performed with the strain 
relief loop prior to this step. Figure 3-20 shows the "N” contact soldering operation with 
the solar cell and interconnector held in the correct relative position in a soldering fixture. 
During process development, joints were made at several soldering iron temperatures in 
order to determine the effect of this parameter on the resulting contact pull strength. 

Table 3-5 summarizes the results of these measurements for a 45 degree puU angle be- 
tween the cell and the interconnector strip. A soldering iron temperature of 310°C (590°F) 
was found to give satisfactory results. In all cases joint failure occurred by a separation 
of the cell metalization from the silicon wafer. 

Following ”N" contact soldering the cells are cleaned untrasonically to remove traces of 
flux. 

Soldering of the "P” contact at three places per cell occurs following bonding of the 19- 
cell array to the glass coverplate. No positioning fixture is reqiaired for this operation 
since the cells are securely bonded to the glass with the PVB film. 
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Figure 3-20. Soldering of Solar Cell ’*N” Contacts 
Table 3-5. 45® Contact Pvill Test Results on Spectrolab Cells 


Solder 

Type 

Soldering 

Iron 

Temperature 

(®F) 

"N” Contact 
Pull Strength 
(grams) 

Remarks 

Sn 60 

700 

397 

20% Ag/Si separation 

Sn 60 

700 

397 

80% Ag/Si separation 

Sn 62 

640. 

1049 

30% Ag/Si separation 

Sn 62 

640 

368 

100% Ag/Si separation 

Sn 62 

590 

1049 

Cell broke 

Sn 62 

590 

822 

Interconnect broke at hole 

Sn 62 

590 

709 

80% Ag/Si separation 

Sn 62 

590 

737 

100% Ag/Si separation 

Sn 62 

590 

567 

100% Ag/Si separation 
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3.3.3 SOLAR CELL BONDING 

The bonding of the solar cells to the embossed surface of the tempered glass coverplate 
is the next assembly step. The 19 cells are loaded in the module bonding fixture shown 
in Figure 3-21. This fixture, which has a pocket for each cell, controls the linear 
positioning of each cell within the modxile as well as the angular orientation of each cell. 

A series of holes drilled into an internal cavity permits the later application of vacuum 
to the space between the fixture and the glass coverplate. A thermocouple has been 
embedded within the fixture to enable monitoring of plate temperature during the sub- 
sequent bonding operation. 

Two disks of 1.75-inch diameter SAFLEX SR-10 polyvinyl butyral are next centered on 
the top of each solar cell. One surface of this product is grooved to enhance the escape 
of air during vacuum lamination. This surface of each of the two disks is placed in con- 
tact with the glass or solar cell. 

A glass coverplate is positioned on top of the cells and PVB disks and taped to the bond- 
ing fixture around the perimeter with high temperature Teflon tape. This fixture is then 
placed on a preheated hot plate as shown in Figure 3- 2 2(a). An insvilation blanket is 
placed on top of the glass to reduce heat loss from the coverplate surface as shown in 
Figure 3-22(b). A thermocouple is placed at the center of the module between the glass 
coverplate and insulation blanket. The insulation blanket reduces the temperature dif- 
ference between the fixture and glass coverplate and produces a more uniform flow-out 
of the PVB from cell-to-cell across the module. 

Table 3-6 gives the temperature- vacuum profile which has been developed to give satis- ' 
factory bonds with two disks of SAFLEX SR-10 as described above. Figure 2-23 shows 
front and rear views of the glass coverplate for Module Serial No. SM-2 following cell 
bonding and soldering of the interconnectors to the "P” contacts. 

Initial attempts to bond cells to the glass coverplate with a single disk of 2 . 00-inch dia- 
meter SAFLEX PT-10 proved to be unsuccessful in terms of the repeatability of obtaining 
a void- free lamination. The thickness of the PVB was doubled by adding another 2.00- 
inch disk of PT-10 to each cell. This approach was successful in eliminating the large area 
voids in the lamination, but bubbles continued to be a problem. Attempts were made to 
dissolve these bubbles by subjecting the coverplate assembli6s to an autoclave cycle which 
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Figure 3-21. Loading of Bonding Fixture 


L 




O ' 


(a) Without Insvilation Blanket 


(b) With Insulation Blanket 
Figure 3-22. Solar Cell Bonding Operation 
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Table 3-6. Typical Cell Bonding Cycle 


Elapsed 

Time 

(minutes) 

Fixture 

Temp. 

(®F) 

Coverplate 

Temp. 

(®F) 

Vacuum 

(in. Hg Gauge) 

Remarks 

0 

75 

75 

0 

Place on hot plate 

6 

150 

110 

5 


14 

200 

170 

10 


16 

220 

190 

20 


22 

250 

225 

27 


38 

285 

266 

0 

Vent for 2 minutes 

40 

285 

268 

27 


84 



0 

Remove from hot plate and cool 


consisted of a pressure of 180 psi at 300®F. The coverplates were maintained at this 
temperature and pressure for at least 30 minutes before the temperature was allowed to 
decrease. The pressure was maintained until the coverplate temperature was below 120°F. 
This cycle was successful in eliminating entrapped bubles, but in some cases these 
bubbles reappeared following overnight soak at ambient conditions. Some yellowing of the 
PVB where it runs out around the circiunference of the cells was also a resvilt of the auto- 
clave cycle. Since this autoclave step is an expensive addition to the processing sequence, 
efforts were made to obtain bubble-free bonds directly. To this end, the change to SAFLEX 
SR- 10, which has more plasticizer and better flow characteristics at a lower temper- 
ature, has resulted in satisfactory bonds without the need for an autoclave pix>cessing 
step. Experiments with the diameter of the SR-10 disks, which included 2.00, 1.87, and 
1.75 inch diameters, has led to the selection of two 1.75-inch diameter disks. This 
selection provides nearly bubble-free bonds with limited run-out around the cell cir- 
cumference and virtually no underrun. 

3.3.4 SUBSTRATE LAMINATION 

The lamination of the semi-flexible substrate, which is the next step in the processing 
sequences , involving bonding the outer FLEXSEAL skins to the epichlorohydrin foam core 
and to the two-sided flexible printed circuit board. The B.F. Goodrich FLEXSEAL con- 
tact adhesive system, which was selected for this application, has been found to give 
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(a) Front 



(b) Back 

Figure 3-23. Coverplate /Solar Cell Assembly for Module Serial No. SM-2 





satisfactory results if care is exercised when the initial contact is made between the ad- 
hesive-coated surfaces. Waxed paper is slipped between these surfaces so that proper 
alignment can be achieved before contact is made. ThiS' paper interlayer is then gradu- 
ally pulled out while the surfaces are rolled together along a line to minimize the amount 
of entrapped air pockets. 

The flexible printed circuit board substrate material has been made large enough to 
cover the entire area of the foam core. This change was found necessary to reduce the 
formation of bubbles luider the FLEXSEAL skins when the mod\ile is exposed to elevated 
temperatures. 

3.3.5 MODULE ENCAPSULATION 

The encapsulation of the module and the bonding of the rear TEXTOLITE cover are both 
accomplished using RTVll. The rear side of the glass coverplate and solar cells are 
primed with SS-4044 followed by a spray coating of RTV108 (10 to 20% by weight mixed 
with heptane). The RTVll is mixed, deaerated and poured over the rear surface of the 
glass coverplate assembly. The material is spread over the entire glass surface until the 
height of the silicone just covers the interconnector strain relief loops. The spread RTV- 
ll is again deaerated in a large vacuum enclosure to remove entrapped bubbles before 
the rear TEXTOLITE cover, which has also been previously primed with SS-4044 and 
RTV108, is attached. 

3.4 MODULE ELECTRICAL PERFORMANCE 

Table 3-7 summarizes the electrical performance measurements made using the Large Area 
Pulsed Solar Simulator (LAPSS) with JPL-supplied Terrestrial Secondary Standard No. 

025 as the reference cell. In the case of the qualification modules which are designated 
by a ”Q” following the serial number, the electrical performance reported in this table 
represents the final measurement made after all environmental exposures. Figure 3-24 
is an example of a typical I-V characteristic as obtained using the LAPSS . The mean 
performance of these modules is 5.79 watts at the maximum power point with a standard 
deviation of 0.169 watts. Module Serial Number SM-7Q was eliminated from the average 
since a cracked cell within this module results in a significant degradation in performance. 

This average modvile performance reflects a net enhancement of 7.8 percent compared to 
the measured average bare cell performance which is 282.6 mW per cell or 5.37 watts 
per module. As a further verification of this resixlt, modvde Serial Number SM-12 was 
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Table 3-7. Svunmary of Modxile Electrical Performance Measurements 


Nodule 

Serial NiMber 

*SC 

(-A) 

*0C 

(Veits) 

(Volts) 

(Watts) 

SN- IQ 

704 

11.4 

9.35 

5.71 

- 20 

714 

11.5 

9.50 

6.16 

- 3 

702 

11.7 

9.30 

5.90 

*Q 

712 

11.6 

9.50 

6.14 

- SA 

710 

11.5 

9.50 

6.05 

- 60 

713 

11.5 

9.50 

5.93 

- 70 

651 


10.10 

4.39 

-10AQ 

695 


9.40 

5.80 

-12A 

713 

11.7 

9.30 

5.99 

-13A 

702 

11.5 

9.50 

5.96 

-17A 

710 

11.5 

9.30 

5.94 

-18 

702 


9.40 

5.88 

-19 

697 

11.6 

9.30 

5.83 

-20A 

697 

11.4 

9.30 

5.90 

-21 

706 


9.40 

6.02 

-220 

706 


9.40 

6.03 

-23 

704 

11.6 

9.30 

5.90 

-24 

692 

11.5 

9.40 

5.85 

-26 

710 

11.5 

9.40 

6^.04 

-27 

714 


9.30 

6.02 

-28 

695 


9.20 

5,77 

-29 

708 

11.6 

9.30 

5.97 

-30 

699 

11.4 

9.25 

5.83 

-31 

696 


9.30 

5.79 

-32 

685 

11.4 

9.40 

5.79 

-34 

699 


9.20 

5.ro 

-35 

703 


9.30 

5.78 

-36 

700 


9.20 \ 

5.72 

-37 

693 


9.30 

\ 5.74 

-38 

698 

11.4 

9.40 

5.70 

-39 

687 

11.4 

9.40 

5.71 

-41 

682 

11.3 

9.35 

5.64 

-42 

703 

11.4 

9.40 

5.84 

-43 

705 


9.40 

5.87 

-44 

684 

11.4 

9.50 

5.86 

-45 

700 

11.4 

9.30 

5.65 

-46 

704 


9.20 

5.72 

-47 

700 

11.4 

9.40 

5.67 

-48 

691 

11.3 

9.20 

5.58 

-49 

687 


9.20 

5.63 

-50 

699 

11.3 

9.25 

5.70 

-51 

698 

11.3 

9.30 

5.67 

-52 

704 


9.25 

5.62 

-53 

695 

11.3 

9.20 

5.56 

-54 

680 

11.0 

9.00 

5.49 

-55 

711 

11.3 

9.20 

5.S6 

-56 

694 


9.20 

5.56 

-57 

697 

11.3 

9.35 

5.62 

-58 

686 

11.3 

9.20 

5.49 

-59 

701 


9.20 

5.5S 


* At 100 iM/ca^ using Term trial Secondary Standard Mo. 02S as tlw reference 
cell. Cell tcaperature corrected to 28*C. 
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Figure 3-24. I-V Characteristic for Module Serial Number SM-4 

built using solar cells which had been previously serialized and illuminated to determine 
the I-V characteristic for each cell. The use of these cells permits the accurate deter- 
mination of the absolute gain due to covering and encapsulation. 

Table 3-8 gives the voltage reading for each cell at four current values which bracket 
the maximum power point of the cells. The sum of the values for all 19 cells at each of 
these four currents yields the composite modvile characteristic with no coverplate. The 
overall module performance with bare cells is 5.56 Watts at a maximum power point 
voltage of 9.25 Volts. This equates to an average bare cell efficiency of 13.3 percent 
based on a nominal total solar cell area of 419.2 cm . 

Following the bonding of these cells to the embossed surface of an ASG SUNADEX glass 
coverplate with PVB film , the I-V characteristic of the entire module was measured be- 
fore encapsxilation with RTVll. Figure 3-25 shows this characteristic both before and 
after encapsulation with the white silicone pottant. Prior to encapsulation the enhance- 
ment in maximum power output is calculated to be 2.3 percent based on the composite 
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Table 3- 8, Performance of Individual Solar Cells in Module Serial No, SM-12* 


Cell 

Position 

CeU 

Identification 

Number 

Voltage at Specified Current (Volts) 

575 mA 

600 mA 

625 mA 

650 mA 

1 

580-1 

0.497 

0.485 

0.466 

0.432 

2 

590-1 

0.504 

0.488 

0.460 

0.350 

3 

590-2 

0.497 

0.485 

0.467 

0.440 

4 

590-3 

0.504 

0.490 

0.468 

0.420 

5 

600-1 

0.508 

0.493 

0.470 

0.410 

6 

600-1 

0.512 

0.500 

0.485 

0.452 

7 

590-8 

0.503 

0.487 

0.460 

0.400 

8 

570-13 

0.495 

0.480 

0.456 

0.400 

9 

600-2 

0.508 

0.495 

0.473 

0.420 

10 

610-1 

0.502 

0.492 

0.474 

0.442 

11 

610-2 

0.507 

0.495 

0.476 

0.437 

12 

580-2 

0.495 

0.480 

0.454 

0.390 

13 

580-5 

0.500 

0.487 

0.460 

0.385 

14 

620-1 

0.499 

0.486 

0.469 

0.435 

15 

580-3 

0.500 

0.488 

0.470 

0.423 

16 

580-7 

0.499 

0.482 

0.452 

0.350 

17 

580-9 

0.499 

0.485 

0.456 

0.375 

18 

580-12 

0.498 

0.484 

0.465 

0.407 

19 

590-10 

0.500 

0.483 

0.455 

0.375 

Total Voltage 

9.52 

9.27 

8.84 

7.74 

Measured on SM-12 





Before Enct^sulation 

9.72 

9.50 

9.00 

8.00 

Measured on SM-12 





After Enc^sulation 

10.00 

9.80 

9.30 

9.22 


* At 100 mW/cm^ insolation and at 28°C 
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Figure 3-25. I-V Characteristics of Module Serial No. SM-12 
Before and After Encapsulation 

characteristic of the 19 bare solar cells. Following encapsulation, this enhancement is 
increased to 7. 7 percent compared to the same composite bare cell characteristic. This 
agrees closely with the result sighted above based on the average performance of all modules. 

It should be noted that this is the net enhancement, which includes the transmission and 
reflection losses associated with the glass and PVB. The exposed surface of the glass cover- 
plate was untreated. The addition of a suitable anti- reflective coating to this outside surface 
could further augment the module output by an additional 1. 5 to 2.0 percent. 

This enhancement in module electrical output can be attributed to the internal reflection of 
light incident on the embossed surface of the SUNADEX glass in the wliite interstices. Ibis 
light, which would be otherwise wasted, is internally reflected within the glass coverplate 
until some portion is absorbed by the solar cell active surface or converted into electrical 
output. This phenomenon which has been called the "zero depth concentrator", was reported 
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under the New Technology provision of this contract (See Section 6). Figure 3-26 shows the 
enhancement in short-circuit current as a function of packing factor. These results were 
obtained using the test specimen shown in Figure 3-27 and varying the size on the circular 
aperture which was placed over the glass coverplate to achieve differing exposed annulus 
areas aroimd the solar cell. Appendix A describes an analytical model for the zero depth 
concentrator enhancement phenomenon. 

3.5 MODULE QUALIFICATION TESTING 

3.5.1 ELECTRICAL PERFORMANCE RESULTS 

The electrical performance of each qualification module was measured initially and again 
after each of the environmental exposures. Table 3-9 summarizes the results of these 
illumination tests which were performed using the LAPSS with JPL-supplied Terrestrial 
Secondary Standard No. 025 used as the reference. 

With the exception of module serial number SM-IOAQ, the change in maximum power output 
due to the thermal cycling exposure was positive with an average improvement of 2 percent. 
This improvement was due to an increase in the curve fill factor (CFF) in all cases, which 
would tend to indicate a decrease in contact series resistance or an increase in shunt resist- 
ance. In the case of module serial number SM-IOAQ the ^pearance of a relatively large 
number of bubbles within the PVB caused a 2. 5 percent decrease in the short-circuit current. 
However, even in this case the module experienced an increase in the CFF as a result of 
the thermal cycling exposure. 

The humidity- temperature exposure was found to cause no significant change in the module 
electrical performance. Similarly the mechanical integrity test did not cause a change in 
module electrical performance beyond the expected accuracy of the measurement. However, 
in the case of module serial number SM-7Q the 22 percent decrease in maximum power out- 
put can be attributed to a single cracked ceU which resulted from a highly localized bearing 
load Imposed by standing on the glass coverplate as discussed in Section 3.5.4. 
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Figure 3-26. Enhance Module Ou4>ut with an Embossed Glass Coverplate 
and White Interstices 



Figure 3-27. 


SUNADEX Glass-Covered Test Specimc 
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Table 3-9. Summary of Qualification Module Electrical Performance 



Initial Measurement 

After Thermal Cycle Test 

After Humidity- Temperature 

1 Test 

After 

> Mechanic 

al InUgrity 

' Test 

Module 

Serial 

isc 

Voc 

VmP 

Pmax 

Isc 

Vex: 

VmP 

pmax 

Isc 

Vex: 

VmP 

Pmax 

Isc 

Vex: 

VmP 

Pmax 

Number 

(mA) 

(Volts) 

(Volts) 

(Watts) 

(mA) 

(Volts) 

(Volts) 

(Watts) 

(mA) 

(VolU) 

(VolU) 

(Watts) 

(mA) 

(VolU) 

(VolU) 

(Watts) 

SM-IQ 

705 

11.3 

9.10 

5.32 

705 

11.5 

9.30 

5. 67 

705 

11.5 

9.30 

5.69 

704 

11.4 

9.35 

5.71 

SM-2Q 

720 

11.7 

9.30 

6.05 

714 

11.6 

9.40 

6.10 

715 

11.65 

9.40 

6.11 

714 

11.5 

9.50 

6. 16 

SM-4Q 

713 

11.6 

9.30 

6.01 

713 

11.7 

9.40 

6. 08 

713 

11.65 

9.35 

6.09 

712 

11.6 

9.50 

6. 14 

,SM-6Q 

719 

11.7 

9.25 

5.79 

715 

11.5 

9.40 

5. 86 

716 

11.6 

9.40 

5.90 

713 

11.5 

9.50 

5.93 

SM-7Q 

707 

11.5 

9.25 

5.68 

70S 

11.5 

9.30 

5.76 

707 

11.6 

9.20 

5.66 

651 

11.4 

10. 10 

4.39 

SM-IOAQ 

713 

11.7 

9.30 

5.92 

695 

11.7 

9.40 

5. 79 

695 

11.6 

9.40 

5.80 

— 

— 

— 


SM-22Q 

706 

11.6 

9.30 

5.91 

712 

11.6 

9.30 

6.03 

712 

11.6 

9.35 

6. 02 

706 

11.5 

9.40 

6. 03 


Notes: (1) Module SM>22Q subjected to 17 thermal cycles. 

(2) Module SM-7Q performance indicates broken cell resulting from bearing load 
due to 68 kg man standing on one leg near center of glass coverplate. This is only 
module which was subjected to this loading. 

(3) Module SM-IOAQ was not subjected to the mechanical integrity test. 




Module serial number SM-IOAQ was not subjected to the mechanical integrity test since 
there were only six positions for qualification modules within the simulated roof test matrix. 

3. 5. 2 THERMAL CYCLE TESTING 

TTie qualification modules were arranged within the thermal cycling chamber as shown in 
Figure 3-28. An aluminum foil tent was constructed over the modules to reduce the rate of 
change of temperature and provide a uniform temperature distribution in the gas surrounding 
the modules. Thermocouples within this tent were monitored during the test and were indica- 
tive of the chamber temperature as the cycling was performed in accordance with the profile 
shown in Figure 3-29. 

A visual examination of the modules, which was performed after the completion of a few 
thermal cycles, revealed that delamination was occurring between the RTVll encapsulant 
and the glass coverplate. This was generally localized around the perimeter of the glass 
coverplate. These modules had a thin coating of diluted RTV108 applied to the rear side of 
the glass coverplate to act as a primer for the RTVll. In view of this development an effort 
was made to improve this adhesion by appl 3 dng SS-4044 primer to the glass coverplate prior 
to coating with diluted RTV108. A module (Serial Number SM-22Q) was fabricated with this 
two step priming system and introduced into the thermal cycling chamber at the conclusion 
of cycle number 33. 

Post-test examination of this module revealed no evidence of delamination resulting from the 
17 cycle exposure while the other six initial modules had shown delaminations after a few 
thermal cycles with no observable increase in area affected after the completion of the 50 
cycles. Based on these observations it can be concluded that the two-step priming system 
provided the necessary adhesion between the glass and the RTVll. 

The formation of bubbles between the top FLEXSEAL skin and the foam core was also ob- 
served at the conclusion of this exposure. The occurrence of these bubbles, which is not 
considered to be a functional problem, appears to be related to the initial uniformity of the 
layer of CA1056 contact adhesive which bonds this interface. A thin, uniform application of 
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Figure 3-28. Test Set-up for Thermal Cycling Exposure 
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TIME (hr) 

Figure 3-29. Thermal Cycle Test Profile 

this adhesive is necessary to prevent bond failure during thermal cycling. Efforts have been 
made on subsequent modules to achieve the necessary uniformity of (his adhesive and to pre- 
vent the initial entrapment of air during lamination. 

3.5.3 TEMPERATURE-HUMIDITY TESTING 

The qualification modules were subjected to tiie temperature-humidity profile shown in 
Figure 3-30. A post-test visual examination of these modules revealed a failure of the bond- 
ed joint between the glass coverplate and the top FLEXSEAL skin. Ibis failure was a lack 
of adhesion between the CA 1056 contact adhesive and the glass surface. These joints were 
reworked by bonding with RTV102 on all modules and using RTV102 as the bonding adhesive 
for all new modules. 
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Figure 3-30. Temperature-Humidity Test Profile 
3.5.4 MECHANICAL INTEGRITY TESTING 

The mechanical integrity testing consisted of 100 cycles of bidirectional loading at 2394 Pa 
(50 lb/ft2) applied as a uniformly distributed load over the exposed coverplate surface of 
each individual qualification shingle. This loading was accomplished as shown in Figure 3-31. 
A vacuum fixture was designed to fit over the top of the glass coverplate. This fixture was 
weighted to obtain a total mass of 12.38 kg (27.3 lb). With an exposed module area of 
0. 0507 m2 this amounts to the required 2394 Pa in the downward direction. The upward 
loading was then obtained by pulling up on the loaded vacuum fixture until the force gauge 
read ^proximately 243N (54.6 lb). This loading cycle was repeated 100 times for each of 
the six qualification shingles. 

The maximum upward deflection measured at the bottom edge of the coverplate ranged from 
25 to 13 mm (1. 0 to 0. 5 inches) depending on the location of the module within the matrix. 
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Figure 3-31. Test Set-up for Mechanical Integrity Testing 
(Photo No. VF 78-249) 
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There was no visual damage as a result of this loading test. However, a post-test examina- 
tion of module no. SM-7Q revealed the presence of one cracked solar cell. Since tins was 
the module which was previously loaded by standing on the glass coverplate as shown in 
Figure 3-32, it is probable that this crack was caused by this highly localized bearing load 
which was at one time 667 N (150 lb) distributed over the area of the ball of one foot. 


3-41 



Figure 3-32. Standing on Installed Shingle Modules 
(Kioto No. VF 78-237A) 
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CONCLUSIONS 


The shingle module design described herein offers many advantages for applications which 
require the mounting of photovoltaic modules on the sloping roofs of new or existing 
residential or commercial buildings since no changes in conventional roof construction are 
imposed by the photovoltaic installation. The ancillary function of this module, viz, , that 
of a weafliertight roof covering; affects additional economy by eliminating the need for a 
conventional roof surface under the photovoltaic module installation. The relatively high 
areal specific output of this shingle module makes its use particularly attractive for installa- 
tion on area limited roof surfaces. An average module maximum power output of 98 watts/m^ 
of exposed module area at the Standard Operating Conditions has been calculated based on the 
measured performance of the 50 modules delivered under this contract. This high specific 
power output can be attributed to the efficient packing of fixe circular ceUs within fixe hexagon 
shape and to the enhancement of the output due to fixe reflected light from the embossed 
glass pattern in the white interstices. An enhancement of 7. 7 percent, compared to bare 
cell performance, has been measured on a shii^e module of Axis design. Another way of 
looking at this enhancement is that the white interstices, which amount to 87 em^ per module, 
perform as solar cells with an efficiency of about 5 percent. 

A typical installation of 1872 modules which are arranged as three circuits, each with 
26 series-connected modules and 24 parallel-connected modules , would require a roof 
length of 10.128 m. A slant height of 9.621 m is required to accommodate three cir- 
cuits in this direction. Thus, fixe required roof area of 97.44 m^ for this typical Installa- 
tion results in an overall area utilization of 0. 805, which is defined as the ratio of solar 
cell area to rectangular roof area. This installation will produce a rated output of 9. 2 kW 
at 205 vdc at the Standard (Iterating Conditions -v^ich include a calculated NOCT of 61°C. 

The installation of a shingle solar cell module roof will require special precautions and pro- 
cedures to protect gainst electrical shock hazards and mechanical damage to the modules. 
Each module could be covered with an opaque pressure-sensitive layer over the exqxosed 
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glass coverplate that would remain in place until after the electrical installation is 
completed, or, as an alternative, the module installation could be performed at night. The 
localized bearing loads associated with standing or kneeling on the installed module cover- 
plate must be distributed over a large coverplate surface area to prevent breakage of solar 
cells due to excessive deflection near the center of the coverplate. Specially constructed 
platforms which distribute the load uniformly over the area of four coverplates should be 
employed during installation to prevent breakage of solar cells. 
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SECTION 5 


RECOMMENDATIONS 


An outdoor exposure test of an assembly of shingle modiiles would seem appropriate as the 
next step.in the application of this development. This exposure is currently underway at the 
JPL Field Test Site, 


Ihe matching of individual solar cells to fully utilize the electrical output of each cell within 
a module is an important consideration which is often overlooked and could have a significant 
impact on the overall cost per watt of delivered module output. This matching requires more 
than just knowledge of the individual solar cell current output at a specified test voltage. Of 
equal importance is the short-circuit current (Isc) of each cell since a cell with low 

IgQj when connected as a series element with cells of hi^er Isc will limit the total module 
output Isc to its value. Therefore, it is recommended that both the Isc aod the current at 
a specified test voltage, which is slightly lower than the anticipated maximum power voltage, 
be measured for each individual solar cell. Matching of cells within a module should be 
based primarily on the Isc readings with secondary emphasis given to the difference between 
the two current readings which is measure of the slope of the "constant" current portion of 
the I-V characteristics. 

Design optimization studies using the analytical model for the zero depth concentrator should 
be performed to determine if further enhancements in module output are possible. 
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SECnON 6 
NEW TECHNOLOGY 


The following New Technology disclosures have been reported to JPL: 


1. Descriptive Title; A Zero Depth Solar Photovoltaic Concentrator 
Date Submitted; November 14, 1977 

Name of Innovator: N. F. Shepard, Jr. 

References : Quarterly Report No. 2-DOE/ JPL 954607-78/2, pages 3-14, 3-15, 
6-1 

2. Descriptive Title; hiterconnector for Overlapping Solar Cell Modules 
Date Submitted; ~March 23, 1978 

Name of Innovator: N. F. Shepard, Jr. 

References; Quarterly Report No. 2-DOE/JPL 954607-78/2, pages 3-11, 3-13, 
3-29, 3-32, 3-33, 3-34 

3. Descriptive Title; Embossed Glass-Covered Shingle Solar Cell Module 
Date Submitted; March 23, 1978 

Name of Innovator; N. F. Shepard, Jr. 

References; Quarterly Report No. 2-DOE/JPL 954607-78/2, pages 1-1, 1-2, 3-1 
through 3-16 
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APPENDIX A 


ANALYTICAL. MODEL OF A ZERO DEPTH 
SOLAR PHOTOVOLTAIC CONCENTRATOR 


A.l INTRODUCTION 

Light incident on the white interstices of an embossed glass covered module results 
in an enhancement of the electrical output. The additional task activity described 
herein concerns the development of an analytical model to calculate the improvement 
possible for a variety of module design parameters. This model is capable of 
assessing the impact on optical performance of the following variables acting inde- 
I pendently : 

1. indices of refraction of each of the encapsuiant layers 

2. cover plate thickness 

3 . cell diameter 

4. cell spacing 

5. surface roughness of outer coverplate surface ' 

6. surface roughness of inner coverplate surface 

7. reflection properties of the interstices 

8. optical transmission of each of the layers through which solar energy 
passes 

9. varying angles of light incidence upon the outer coverplate surface 

This model will permit the optical performance of a glass-covered, hexagon-shaped shingle 
solar cell module, in terms of the enhancement in the annual energy incident on the 
active solar cell surfaces, to be evaluated for various combinations of the above mentioned 
variables . 

In combination with algorithms which define the shingle module cost as a function of 
cell diameter and spacing, it will be possible to evaluate the shingle module cost per 
unit of annual energy output as a function of shingle geometry. 
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A. 2 GENERAL DESCRIPTION OF OPTICAL MODEL 

The factors to be considered in the optical model can be visualized with the aid of Figure 
A-1 which shows a section through the glass coverplate with the underlying solar eeUs. 
Related to the letter designations on the figure, the factors to be considered include the 

following: 

A. orientation of the earth-sun line relative to the normal to the cover- 
plate of the shingle module 

B. surface coating at the air-to-glass interface 

C. structure of the outside or exposed surface of the glass coverplate. 

In the case of SUNDADEX glass, this surface has a slightly matte 
finish. 

D. the glass material 

E . structure of the inside surface of the glass coverplate . In the case of 
SUNADEX glass, this surface has a deeply embossed regular pattern of 
indentations. 

F. the material in contact with the inside surface of the glass coverplate in the 
interstices between the solar cell pattern 

G . the transparent solar cell bonding material 

H. the solar cell spacing and configuration 

As described above, the optical portion of this model entails the consideration of many 
factors which influence the ultimate conversion of incident solar energy into electrical 
energy at the solec? cell array. For maximum annual electrical energy output from a 
shingle, the goal is to maximize the’ useful solar energy which crosses the glass cover- 
plate outside surface and to maximize the containment and absorption, by the solar cells, 
of the admitted energy. 

The analytical approach to the optical portion of this problem will first consider the 
incident solar energy particularly as it relates to the angle of incidence effects. Light 
incident on the shingle module and traversing the glass is definable in terms of two 
angles as shown in Figure A- 2: an angle of incidence (I) , measured with respect to the 
shingle normal; and an azimuth angle (A) , measured with respect to a due south axis in 
the plane of the module surface. 




Figure A-1. Section Through Glass Coverplate Showing Factors 
Influencing Enhancement of Output Due to Zero Depth 
Concentrator Phenomena 
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Figure A-2. Earth-Sun Line Orientation 
on Shingle Module Surface 



Since the solar cell pattern is not point-symmetrical, performance will, in general, 
depend on azimuth, which will vary with time. Over a period of a year^ however, the 
dependence on azimuth wUl be at least suppressed, if not washed out completely. There- 
fore, since the ultimate concern is with energy collected over a one year period, the 
problem will be defined only in terms of incidence angle, with all dependencies on azimuth 
being averaged out in the computation. 

The basic form of the input data.wiU be a histogram representing the annual energy 
incident on the glass coverplate surface over 2 degree intervals of the angle of. incidence . 
As shown in Figure A- 3, this histogram can be represented in the data structure as an 
array W (I) of 45 elements. ■ 



ANGLE OF INCiDENCE INTERVAL DEGREES) 


Figure A- 3. Histogram of Annual Input Direct Beam 

Solar Energy 

The various refraction and reflection processes within, the transparent coverplate are 
shown in Figure A- 4, Process represents the entrance of incident solar energy 
into the coverplate. A portion of this energy is lost by a Fresnel reflection process 
at the rare to dense outer surface interface. Figure A- 5 represents a typical Fresnel 
reflection process function. 





This function can be modified by appropriate anti-reflective surface coatings. The 
specular beam of energy which passes into the coverplate is refracted ^by a Snell s 
Law refraction process where the ray is deflected toward the surface normal when pass- 
ing from the rare-to-dense medium. The exit angle, 0 is given by; 


0 = sin 


- 1 /sin <l> 




Similar reflection and refraction processes occur at the internal coverplate-to-solar 
cell interface , represented by Process ; at the internal coverplate -to-interstices 

interface, represented by Process (|) ; and at the Internal coverplate-to-air inter- 
face, represented .by Process ® • 


At each of these interfaces the sp.ecular, beam is modifed by a scattering function which 
is shown diagramically in Figure A“^6.- 



Figure A-6. Scattering at'Outer and Inner 
Coverplate Surfaces 


A scattering function , which operates on W ( 0 ) by convolution , can be identified for 
each of these surfaces to account for this angular spread of the incident beam. 



The passage of light through the coverplate material is also influenced by a bulk 
absorption process which is shown in Figure A-7. The absorption coefficient k is 
defined as an average over the sunlight spectrum. 

Beam division processes also play an important role in the analytical definition of optical 
performance of the zero depth concentrator. Figure A- 8 shows an array of circular 
beams representing reflected light from a solar cell pattern at an inclination angle <i». 

In making two passes through the coverplate (bottom to top and back to bottom) this 
pattern is sheared with respect to the original by an amount , S , given by : 

S = 2t . tan <t> 

The overlap of the sheared pattern and original pattern is defined as the beam division 
function AA(s). Compound beam division functions are definable as a consequence of 
sequential passes through the coverplate as shown in Figure A- 9. 

A diagram of the various processes involved in the optical model is given in Figure A- 10. 
Starting at the top of the page, this diagram shows the sequential process steps which 
operate on the incident energy distribution array W (I) to produce consecutive definitions 
of W as the energy passes through the coverplate for the first time and is incident on a 
solar cell, in the case of the ”A" branch of the flow diagram, or on the interstices, is 
the case of the "A" branch. The various processes described above and enumerated on 
Figure A-10 operate on each "2-pass” branch after the initial reflection and beam division 
at the bottom surface of the coverplate. The "2-pass" sequence of operations relates 
to energy which passes from the bottom surface through the glass , is reflected and 
scattered from the top surface, and passes through the glass again before being incident 
on the bottom surface. At each occurence of reflection at the bottom surface an areal 
beam division, takes place , accompanied by an associated translation or shear of the geo- 
metric pattern of circles within a hexagon. The program continues until the light 
energy has reached the bottom surface for the third time, or until two "2-pass" sequences 
have been completed. At appropriate places within the sequence the amount of energy 
incident on the solar cells is accumulated. The amount of energy available for collection 
is successively reduced as the chain of processes continues until after two "2-pass" 
sequences it is felt that the residual energy will be quite small. 
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AA(Si) A(S2) 




Figure A-9. Definition of Beam Division Functions 
A. 3 ASSUMPTIONS 

The development of the analytical model will be based on the existing shingle solar ceU 
module concept which consists of an array of circular solar cells arranged in a hexagon 
configuration under. a common glass coverplate, as shown in Figure A-11. The cell 
diameter (d), spacing (S) and the number of solar cells (7 or 19) will be varied in the 
model but the basic hexagon geometry will be retained. For this module geometry the 
edge dimension of the hexagon (E) can be expressed in terms of the aforementioned 
variables by the following relationships. 

E = 2.57735d + 3.1547S (for 19 cell module) 

E = 1.57735d + 2.1547S (for 7 cell module) 

Using the shingle module geometry shown in Figure A-12, the areas of the coverplate. 

(A ) and substrate (A ) can be expressed as a function of E by the following ex- 

s 

pressions : 

A ^ = 2.598E^ 
cp 

A = 3.8971E^ 
s 
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Figure A- 10. Diagram of Processes Involved in the Optical Model 
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The minimum required glass coverplate thickness is given in Figure A-13. as a function 
of the hexagon edge dimension. A constant deflection at the center of a simply supported 
plate subjected to a central point load was used as the criterion for establishing this 
curve. It was assumed that the existing shingle design) which uses 3.2mm thick glass 
with a 139.7 mm edge dimension, gives a marginally acceptable deflection under these 
loading conditions. 

Only direct beam radiation will be considered in the model. The assumption is made that 
the combination of factors which optimize the enhancement due to the direct component 
of insolation will also act in a similar manner to enhance the contribution due to the 
diffuse component. Since the angle of incidence on the outer coverplate surface has a 
major influence on the enhancement due to this phenomenon, the evaluation of the system 
performance will be made based on calculated hourly values for the angle of incidence 
and corresponding direct beam insolation component on the module surface. SOLMET 
■data tapes for Boston, MA and Phoenix, AZ will be used to represent two diverse sites 
for this analysis. Figures A- 14 and A- 15 give the accumulative distribution of the 
annual direct beam energy density on the module surface as a function of angle of 
incidence for Bostom, MA and Phoenix, AZ, respectively. Hourly SOLMET data were 
used in both cases with a roof slope angle of 10 degrees less than the site latitude. 

These data show that a fixed E-W tilted solar array collects more than 50 percent of its 
annual direct energy at incidence angles of greater than 30 degrees. Table A~1 tabulates 
these data in 2 degree angular intervals for use as an input to the model. 

A. 4 EXPERIMENTAL RESULTS 

A. 4.1 SCATTERING FUNCTIONS FOR SUNADEX GLASS 

The determination of the scattering functions for each surface of SUNADEX glass involved 
the measurement of the reflected energy which passes through an aperture of known size. 
Two different experimental test set-ups were required as shown in Figures A-16 and A-17. 
The first of these was used to determine the scattering function for the matte surface 
of the .coverplate, but a focusing lens could not be obtained with a small enough f/ 
number to permit the measurement of the more widely dispersed pattern from the 
embossed surface. Consequently the test set-up shown in Figure A-18 was used for 
this case. 
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Figure A-13. Minimum Kequired Glass Coverpiate 
Thickness vs. Hexagon Size 
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Figure A- 14. Accumulative Distribution of Annual 
Direct Beam Energy Density as a Function of 
Angle of Incidence for a Boston, MA Site Location 


Figure A-15. -Accumulative Distribution of -Annual 
Direct Beam Energy Density as a Function of 
Angle of Incidence for a Phoenix, AZ Site Location 
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Table A-1. Angular Distribution of Annual Direct .Bearo< Energy ^Density for 

Boston, MA and Phoenix, AZ 




BOSTON, 

MA 1960 SOLMET DATA TAPE 

PH0ENI2, A2 

1953 SOLMET DATA TAPE 

ANGLE 

of 

INCIDENCE 

ANNUAL DIRECT 
BEAM ENERGY DENSITY 
ON MODULE SURFACE 
AT ANGLES OF 
INCIOENCE.>0, 

ANNUAL DIRECT 
BEAM ENERGY DENSITY 
ON MODULE SURFACE 
IN INTERVAL 
(0^ - 0?) 

(kWh/m^) 

FRACTION OF ANNUAL 
DIRECT BEAM ENERGY 
DENSITY ON MODULE 
SURFACE IN INTERNAL 
(01 - Ug) 

AHHWAl DIRECT 
BEAM ENERGY DENSITY 
OS MOeUtE SURFACE 
AT ANGLES OF 
INCIDENCE> 0, 
(kWK/m^) ' 

ANNUAL DIRECT . , 
BEAM ENERGY DENSITY 
ON MODULE SURFACE 
IN INTERVAL 

FRACTION OF ANNUAL 
DIRECT BEAM ENERGY 
DENSITY ON H[Mm.E 
SURFACE IN INTERVAL 
(01 - »z) 

Pi 

»2 



(k'Wh/m*^) 

0 

2 

932.4 

0. 

0. 

1853.6 

16.7* 

0.00901 

2 

4 

932,4 

.9*6 

0.01030 

1836.9 

2 o,r. 

0,01084 

4 

6 

922,8 

15.5 

0,01$62 

1816.8 

19.3 

0.01041 

6 

8 

907.3 

, 12.1* 

0.01298 

1797.5 

19.3 

0.01041 

8 

10 

895.2 

12.8 

0,01373 

'1778.2 

ao.s 

0.01106 

10 

12 

882.4 

' 33.9 

0.03636 

1757,7 « 

26.2 

0,01413 

12 

14 

848.5 

40,8 

0.04376 

1731.5 

43.1 

Q.C2325 

14 

16 

807.7 

34,9 

0,03743 

1608.4 

130.2 

0.C7024 

le 

18 

772.8 

37.9 

0.04065 

1558.2 

90.9 

0.04904 

18 

20 

734.9 

46.1 

0.04944 

1467.3 

101.1 

0,05454 

20 

22 

688.8 

- 30.4 

0,03260 

1366,2 

54.1 

0,02919 

22 

24 

658.4 

40.0 

0.04290 

1312,1 

29.4 

0,01586 

24 

26 

618.4 

32,1 

0.03443 

1202.7 

27.4 

0.01478 

26 

28 

586.3 

47,2 

0,05062 

1255.3 

30.0 

0.01618 

28 

30 

539.1 

40.8 

0,04376 

1225.4 

144.8 

0.07812 

30 

32 

498.3 

22.0 

0,02360 

1080,6 

173.0 

0,09333 

32 

34 

476,3 

70.0 

0.07SQ8 

907.6 

101.9 

0.0549? 

34 

35 

406.3 

59.7 

0,06403 

805.7 

68.0 

0.03669 

36 

38 

346.6 

22.9 

0,02456 

737.7 

3S.0 

0.01888 

38 

40 

323.7 

22.2 

0.02381 

702.7 

21.6 

0.01165 

40 

42 

301.5 

54.8 

0,05877 

681.1 

24.2 

0.01306 

42 

44 

246.7 

26,0 

0.02789 

656.9 

106.2 

0,05729 

44 

46 

220.7 

10.6 

0.01137 

550.7 

176,7 

0.C9533 

46 

48 

210.1 

35.8 

0,03840 

374.0 

24.6 

0.01327 

48 

50 

174.3 

34.8 

0,03732 

349.4 

21.4 

0.01155 

50 

52 

139.5 

9.2 

0,00987 

328,0 

17,9 

0.00966 

52 

54 

130.3 

14.3 

0.01534 

‘ 310.1 

20,1 

0,01084 

54 

56 

116.0 

* 28.3 

0.03035 

290. D 

26.5 

O.C1430 

56 

58 

87.7 

12.9 

0.01384 

263.5 

38.1 

0.02055 

58 

60 

74.8 

5.0 

0.00536 

225.4 

103.9 

0.05605 

60 

62 

69.8 

9.5 

0.01019 

121.5 

19.3 

0.01041 

62 

64 

60.3 

24.4 

0,02617 

102,2 

11,4 
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Figure A- 16. Schematic of Experimental Apparatus for the Determination 
of the Scattering Function of the Matte Surface 
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Figure A-17. Schematic of Experimental Apparatus for the Determination 
of the Scattering Function of the Embossed Surface 
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Figure A-18. Experimental Set-up for the Determination of 
Scattering Function of Matte Surface (VF78-526A) 




For the limited scattering from the matte surface the apparatus shown in Figure A~18 
proved to give satisfactory results. A lamp, whose brightness is controlled by a Variac, 
is focused on the pin hole aperture. Energy passing through the pin hole is collimated 
and passes through an iris whose function is to limit the collimated bundle to a controlled 
size. 

Approximately 50 percent of the ihuminating bundle is reflected by the beam splitter 
face. The SUNADEX glass has previously been aluminized with a thin but opaque coat- 
ing. The coating was deliberately kept thin in order to retain' the print of the rough 
surface. The reflected beam, which is diverging due to the diffuse reflection, passes 
through the beam splitter. It is important to keep the bundle as large as possible, Isut 
relatively small with respect to the size of the beam splitter, to prevent internal re- 
flection from interfering with the measurements. 

An "image” of the diffuse bundle is formed at the focal plane of the focusing lens. 

•This plane coincides with the variable size analyzing iris. The diameter of the iris is 
related to the slope by ; 

D = 4f tana 

where 

D = the image diameter 
f = the focal length of the focusing lens 
= 3.00 inches for the test set-up in Figure A-t6 
the surface slope error 

The signal received is a function of the analyzing iris diameter. If the pin hole were 
an infinitesimal point source, and the sample perfectly smooth, the output signal 
would be constant as the analyzing aperture diameter is varied. Since the hole is 
finite, however, there is a variable signal as a function of the analyzing aperture but: 
only for very small values of diameter. If the scattered "image" is fairly large compared 
to the pin hole image, the analyzing function is essentially constant representing a delta 
function input . In such a case the measured signal from the matte surface , as a function 


A- 19 



of iris diameter, is a direct function of surface slope error. The function of the field 
lens is to relay .the energy passing through the analyzing iris to the detector. Its focal 
length is chosen so that the entire diameter of the focusing lens is imaged on the detector 
surface and its diameter is chosen to match the largest surface slope deviation to be con- 
sidered. In this way, all energy passing through the focusing lens is imaged on the 
detector, 

A photovoltaic detector coupled to a lock-in voltmeter whose reference signal was derived 

' , 

from a mechanical chopper was employed-. A chopping frequency of 200 Hertz was em- 
ployed to avoid interference from room light which operate on a base frequency of 120 
Hertz. 

In the case of the embossed surface, shown in the scanning electron micrograph of 
Figure A-19 the test set-up was changed to accomodate the wider scatter of the reflected 
light from this. surface. As shown schematically in Figure A- 17 a laser beam impinges on 
the .embossed glass surface, and is reflected onto a Mangin mirror through an analyzing 
iris. The Mangin mirror then focuses the scattered radiation, which passes through the 
iris, onto a photovoltaic detector which yields an electrical signal proportional to the 
intensity of the light on the detector. The formula for the determination of surface 
slope error is the same as previously presented for the matte surface set-up except 
that (f) becomes the distance between the plane of the analyzing iris and the pdiht of 
focus on the embossed surface. In the case of the set-up. shown in Figure A-17, the 
value of (f) is 4.875 inches. 

Figure A-20 gives the normalzed detector signal as a functioh of analyzing iris diameter 
for both surfaces of the SUNADEX glass. These data -were .fitted with a polynomial of 
the form: 


V = + a^ D agD^ + 

+ agD® + agD® + • 

where 

V = normalized detector .signal 

D = analyzing iris diameter ' (inches) 



NORMALIZED DETECTOR SIGNAL (V) 



Figure A-19. SEM of the Embossed Surface of SUNADEX Glass 



0 25 so 75 100 (MM) 

ANALYSING IRIS DIAMETER (D) 


Figure A- 20. Normalized Detector Signal vs. Analyzing Iris Diameter 



For these two surfaces the values for the polynomial coefficients are .as given in 
Table A-2. 


Table A-2. Coefficients of the Pol 3 momial 
Describing the Detector Signal as a Function of the 
Analyzing Iris Diameter 



Embossed 

Surface 

Matte 

Surface 

^0 

1.44754E-4 

i ^ * 

-4.10052E-3 


1.31156E-1 

6.65302E-1 


-6,=97.379E-1 

1.51415E1 

2 



^3 

■1.47873 

-4.40255E1 

^4 

t1.03T39 ' 

\ 4..56052E1 

^5 

3.44919E-1 

^ . ? 

-1.63851E1 

^6 

-5.,70684E-2 

0 

^7 

3.75892E-3 

0 


The derivative of these expressions yields the fraction of the scattered energy which 
occurs within an infinitesimally small surface slope error. Figure A-21 presents these 
results expressed as a fraction of the surface area which can be characterized by a 
given surface slope error. It should be noted for the case of the matte surface the 
data pertains to a normally incident specular beam while for the embossed surface the 
angle of incidence of the specular beam was approximately 20 degrees.. There is also 
some uncertainty regarding the shape of these curves for surface slope errors of less 
than 10 milhradians due to limitations on minimum opening for the analvzing iris. 

A. 4. 2 BEAM 'DIVISION FUNCTIONS FOR CIRCLES WITHIN A HEXAGON 
The beam division functions for seven circles packed within a hexagon were determined 
experimentally using the test 'set-up shown in Figure A-22. Glass photographic plates 
were prepared for the four different module geometries shown in Figure A- 23, which 
represent overall packing factors (PF) ranging from 0.814 to 0.514. These glass plates 
representing the various patterns were shifted in the test’ setup to obtain the various 
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Figure A-21. Characteristic Reflective Scattering Functions 

for SUNADEX Glass 
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Figure A- 22. Schematic of Experimental Apparatus for the 
Determination of the Beam Division Functions 
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(c) PF = 0.658 (d) PF = 0.514 

Figure A- 23. Photographic Masks used for the Determination 
of Beam Division Functions 






beam division functions given in Figure A-24. The nomenclature DD , DD refers to 
energy which is reflected from a solar cell (D) or from the interstices (D) and con- 
sequently falls on a solar cell (D) after a single reflection from the top coverplate 
surface. This nomenclature is. further- clarified by referring to Figure A-25 which 
graphically rejpresents the four possible combinations with a single reflection from the 
top coverplate -surface . 

For the case where two reflections are involved the nomenclature consist of three 
characters, e.g. , DDD. In this example the beam division function represents that- 
energy which falls on a solar cell after two reflections from the top coverplate surface 
of energy which falls on a solar cell. 

A. 5 DESCRIPTION OF OPTICAL ANALYSIS PROGKAM 

A . 5 . 1 OVERALL PROGRAM FLOVv 

The overall organization of the optical analysis program is shown in Figure A-26. The 
program is written in Honeywell Level 66/6000 FORTRAN and should be compatible, with 
minor modifications, with any FORTRAN compiler. The program is configured for exe- 
cution in a time-sharing mode from a remote terminal. An introductory program explan- 
ation will be outputted on request after compilation has been completed. Input data 
will then be requested before the main portion of the program is entered. As shown in 
Figure A- 27 this, main program cpnsists of an accounting of various integer identifiera 
or flags and repeate'd , calls pf two principal calculation procedures identified as "1-PASS" 
and ■" 2-PASS V.’ ;!.These two prbcedures consist of the sequential call of five subroutines 
as shown in Figures A-28 and A- 29, respectively. 

The details of the main program logic are given in Figure A- 27 along with the status-of 
the five integer, identifiers or flags at each stage of the calculation sequence-. The 
notations D, b„. (D)D ,‘ etc. can be explained with the. aid of Figure A- 30 which 

shows the successive routes, through t,he glass coverplate for the two 2-PASS sequences! 
considered in the model. 

The notation "C" beside some blocks on Figure A- 27 indicates that energy is collected 
by the solar cell as a consequence of that set of processes. 
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\ Figure A- 24. Beam Division Functions 
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Figure A- 25. Graphical Representation of Beam Division Functions 
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Figure A- 26. Overall Program Organization 
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Figure A-28. Logic Flow Diagram for the 1-PASS Procedure 
A. 5.2 SUBROUTINE REFLECT 

Figure A-31 gives a detailed logic flow diagram for the REFLECT subroutine. The 
function of this subroutine is the account for the Fresnel reflection from both the upper 
and lower surfaces of the glass coverplate. The logic diagram, has been annotated to 
indicate the nature of the operation being performed in each block. The routine includes 
the flexibility to consider a single-layer anti-refiective coating on the upper surface of 
the coverplate. The following definitions are appropriate to the understanding of this 
routine: 

F=1 Upper surface, rare-to-dense interface 

F=2 Lower surface, solar ceE interface 

F=3 Lower surface, interstices interface 

F=4 Upper surface, dense-to-rare interface 

N(l) = Index of refraction for air (=1) 


A-33 






















J PD)D 



Figure A-31. Logic Flow Diagram for Subroutine REFLECT 
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Fi^re A-31. Logic Flow -raagram for Subroutine REFLECT (Cont‘d) 
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Figure A- 31; Logic Flow Diagram for Subroutme REFLECT (Cont^d) 






RP:=(COS(l)-(N/M)*SQRT(R))/ 
4(C0S(1) + (N/M)*SQRT(.R)) „ 
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((N/M*COS(l) + SQRT(R)) 
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T:=1-(RP**2+RS**2)/2 
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Figure A- 31. Logic Flow Diagram for Subroutine REFLECT (Cont'd) 
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Figure A-31. Logic Flow Diagram for Subroutine REFLECT (Cont'd) 













Figure A- 31. Logic Flow Diagram for Subroutine REFLECT (Cont'd) 
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N(2) = Index of refraction for glass coverplate 
N(3) = Index of refraction for cell bonding adhesive 
N(4) = Index of refraction for encapsulant 

The routine treats reflection at each of. the interface, identified by by the appropriate 
assignment of the indices of refraction. 

This routine also includes the logic for the accumulation of the; collected energy on the 
solar ceU. 

A. 5. 3 SUBROUTINE REFRACT 

Snell's Law refraction wMeh occurs at the air-to-glass interface as the energy enters the 
glass (Interface F=l) is accounted for by this s.ubroutine as shown in Figure A-32. The 
refraction process changes the angular boundaries of each two degree energy cell. A 
major portion of this procedure deals with the technique employed to redistribute the 
refracted energy into the originally assigned 45 angular energy cells. The refracted 
boundaries (A and C) for each cell are calculated and the redistribution of the energy 
into the original cell boundaries is carried out as shown in Figure A- 33. 

A. 5. 4 SUBROUTINE SCATTER 

The SCATTER subroutine shown in Figure A-34 performs a piecewise convolution of 
the energy function (W) with the scattering function (S) . This latter function takes 
the form of a set of data tables as defined in Figure A- 35. For the SUNADEX glass 
surfaces under consideration the values for these, data tables can be derived from the 
experimentally determined function shown in Figure ‘A- 21. Tables A- 3 through A- 5 
give the required scattering function data tables for the four interfaces in question. 
Both of the embossed (or lower) interfaces will' have the same scattering function table, 
but the matte (or upper) surface interfaces (F=l and F=4) will be different because 
scattering is by refraction in one case and' by reflection in the other. 

A. 5. 5 SUBROUTINE ABSORB 

The ABSORB subroutine shown in Figure A-36 accounts for the broad-band bulk 
absorption which occurs as the energy passes through the glass cover plate material. 
Table A— 6 gives typical values for the broad— band bulk absorption coefficient, G, for 
several soda-lime glass compositions. 
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Figure A" 32. Logic Flow Diagram 


> 


SKIP COATING LAYERS 


INITIALIZE ANGLE COUNTER 


LOAD LOCAL ARRAY C(l) 


> 


RE-INITIALIZE 


M (1-1)1 REFRACT BOUNDARIES 

(1 + 1)1 OF ENERGY CELL 


Subroutine REFRACT 










Figure A- 32, Logic Flow Diagram for Subroutine REFRACT (Cont’d) 
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Figure A- 34. Logic Flow Diagram for Subroutine SCATTER 








Figure A-34. Logic Flow Diagram for Subroutine SCATTER (Cont'd) 
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Figure A-34. Logic Flow Diagram for Subroutine SCATTER (Cont'd) 
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Figure A-35. Definition of Scatter Function Table 


Table A- 3. Scatter Function Table for SUNADEX Glass ~ 

Matte Surface (F=l) 


1 

1 0 

3 0 

(Donor 

Cell) 5 0 

7 0 

9 0 

11 
13 


J (Receiver Cell) 

3 5 7 9 

0 0 0 0 

0 0 0 0 . 

0 0 0 0 

0 0 0 0 

0 0 0 0 


11 13 15 17 

0 

0 

0 

0 

0 
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Table A-4. Scatter Function Table for SUNADEX Glass - 
Matte Surface (F=4) 

J (Receiver Cell) 

i 3 5 7 9 11 13 .15 17 


I 

(Donor 

1 

“.08 

.04 

0 

0 

0 

3 

.04 

-.08 

.04 

0 

0 

Cell) 

5 

0 

.04 

-.08 

.04 

0 


7 

0 

0 

.04 

-.08 

.04 0 


9 

11 

13 




# • • 


19 • * • 
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Table A-5. Scatter Function Table for SUNADEX Glass - 


Embossed Surface (F=2 and F=3) 
J (Receiver Cell) 



1 - 

3 

5 

7 

9 

11 

13 

15 

17 

19 

21 

23 

1 

-.94 

.31 

.34 

.19 

.09 

,01 

0 

0 

0 ■ 

0 

0 


3 

.155 

-;94 

.155 

.34 

.19 

.09 

.01 

0 

0 

0 

0 


I 5 

.17 

.155 

-.94 

.155 

.17 

.19 

.09 

.01 

0 

0 

0 


7 

.095 

.17 

.155 

-.94 

.155 

.17 

.095 

.09 

.01 

0 

0 


(Donor 9 

.045 

.095 

.17 

.155 

-.94 

.155 

.17 

.095 

.045 

.01 

0 


Cell) 

.005 

.045 

.095 

.17 

.155 

-.94 

.155 

• .17 

.095 

.045 

.005 

0 

13 ■ 

0 

.005 

-.045. 

.095 

.17' 

.155 

-..94 

.155 

.17 

.095 

.045 

.005 


15 

17 


« 

9 

9 



Figure A- 36, 


Logic Flow Diagram for Subroutine ABSORB 
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Table A-6. Broad-band Bulk Absorption Coefficients for Several 
Soda-Lime Glass Compositions 


Glass Identification 

Iron Oxide Content (%) 

G (mm ^) 

ASG SUNADEX 

0.01 

0.00184 


ASG LO-IRON 

0.05 

0.00926 


AGS STARLUX 

0.12 

0.01933 



A. 5. 6 SUBROUTINE DIVIDE 

Figure A- 37 gives the logic flow diagram for the DIVIDE subroutine which accounts for 
the beam division which occurs at each reflection from the lower surface of the cover 
plate. An empirically determined beam division function (F) is used to describe the 
geometric fraction of reflected area which is coincident with the solar cells. The deter- 
mination of these functional relationships is described in Section A. 4. 2. 

A. 6 DESCRIPTION OF COST MODELS 

The expected costs- for shingle solar cell modules have been developed as a function of 
solar cell area, cover plate area and substrate area based on the following ground rtiles 


• 500 kW peak 

• Late 1980 time period 

• 13 percent bare ceU efficiency at 28®C 

• Circular solar cells packed within a hexagon shape as shown in Figure A-11. 
Both a 19 cell and a 7 cell configuration will be considered in the cost model. 

• The shingle module configuration is as shown in Figure A- 12. 


The shingle module cost has been divided into three categories as described below: 


1. Solar Cell Cost - The specific cost of the circular solar cells is given in 
Figure A- 38 as a function of the total solar ceU area per module and the 
module configuration (7 or 19 cells). This cost element includes the inter- 
connection between cells as well as the cost of laminating the cells to the 
glass coverplate.' 

2. Encapsulation Cost - The specific cost of the encapsulation system, which 
includes the glass coverplate , the hard back and the solar cell encapsula- 
tion material, is given in Figure A- 3 9 as a function of the coverplate area 
per module and the glass coverplate thickness. 
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Figure A- 37. Logic Flow Diagram for Subroutine DIVIDE 






















SPECIFIC SOLAR CELL COST - INCLUDING INTERCONNECTION AND 






3. Substrate Cost - The specific cost of the semi-flexible shingle module sub- 
strate , which includes the flexible circuit board , the outer sMns and the 
core material, is given in Figure A- 40 as a function of the substrate area 
per module and the glass coverplate thickness. 

Applying the cost model to the existing shingle module geometry with a cell diameter of 
53 mm, a hexagon edge dimension of 139.7 mm, a cover plate area of 0.0507 m and a 
substrate area of 0.0761 m^ 3 rields the results listed in Table A- 7. 


Table A-7. Cost Prediction Using Existing Shingle Mod\jle Geometry 


Item 

Cost per 
Module 
(1975$) 

Fraction of 
Total Cost- 

e Solar Cells (including 
interconnection and 
lamination) 

19.53 

0.834 

e Encapsulation 

2.03 

0.087 

« Substrate 

1.86 

0.079 

Total Cost 

23.42 

1.000 


Such a module* will produce 5-88 watts using a 13.0 percent bare cell efficiency at 28°C. 
At the Standard Operating Conditions, which include a'Calculated NOCT of 61°G‘, this 
module output vWill be reduced to 5.03 watts. Thus, the resulting cost^to-power ratio 
is $4. 66 /peak watt. 

■ A-.7; SAMPLE CASE .CALCULATIONS 

A sample ciculation using the optical analysis program for the zero depth concentrator 

is shown in Table A- 8. This program has been written for time-share execution from a 

remote terminal. In the example given the responses required by the user have been 

underlined. The. results of the calculation for the Phoenix distribution- of incident 

energy as a function of angle of incidence (from Table A-1) are given as the energy 

collected on the solar cells in each two degree angular increment representing the angle 

of the refracted specular ray as it passes through the top coverplate surface. The total 

2 

collected energy is calculated to-be 1644.52 kWh/m . With an input insolation of 1853.6 
kW/m (from Table A-1) this result represents an overall enhancement of 16 percent for 
the packing factor used. 
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Figure A- 40, Specific Substrate Cost 


Table A-8. Sample Calculation Using Optical 
Model of Zero Depth Concentrator 

NO. OF COftTIHG LRYER> 1 OR 0 

INDEXES OF REFPfiCTIDH OF GLRSS COVER » BONDING ADHESIVE %c ENCAtSULANT 
NODULE ^CKI^G' FACTOR ^BETWEEN .512 & .S14> 

GLR^ THICKNESS <:3.£» 4 OR 4.8MM>» BULK ABSORP. COEFF.<l/NM> 

=±> ■ (nns4_ 

REFLECTANCE OF ENCR&SULANT 


0- 2 
19.76 

roro 

1 

■ 

4- 6 
22.22 

6- 8 
29.-33 

8-10 

91.70 

10-12 

117.59 

12-14 

92.47 

14-16 
40. 08 

16-18 

38.19 

18-20 
205. 19 

20-22 

149.28 

22-24 

68.26 

£4-26 

37.14 

£6-28 

£19.25 

£8-30 
76. 16 

3 0-32 
41 . 66 

32-34 

85.89 

34-36 
135. 06 

36—38 

37.63 

38-40 

98.11 

40-42 

14.85 

42-44 

1.27 

44-46 

0.81 

46-48 

0.47 

48-50 
0 . £3 

50-52 
0. 11 

52-54 
0. 04 

54-56 
0. 01 

56-58 
0. 00 

58-60 
0. 00 

60-62 
0. 00 

62-64 
0. 00 

64-66 
0. 00 

66-68 

0. 

68-70 

0. 

70-72 

0. 

72-74 

0. 

74-76 

0; 

76-78 

0. 

78-80 

0. 

80-82 

0. 

82-84 
0- . 

84-86 

0. 

86-88 

0. 

.88-90 

0. 


TOTAL ENERGY COLLECTED= 1644.52 
NAKING MORE RUNS?1=YES» 0=NO 
= 0 
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The model was used to calculate the enchancement', relative to the cosine relationship, 
as a function of angle of incidence on the outer coverplate surface for these same con- 
ditions of module packing factor, etc. The result shown in Figure A- 41 shows a signif- 
icant enhancement in the performance of this type of module relative to the cosine re- 
lationship which is representative of the output of a closely-packed planar module. 

A. 8 FORTRAN CODE 

Table A- 9 is a listing, of the FORTRAN code of the optical analysis program developed 
for the zero-depth concentrator. This pmgram is compatible with the Honeywell L66 
system- but can be executed on any system with a FORTRAN compiler with minor modifi- 
cations to the source statements. 



Figure A- 41. Output Enhancement vs. Angle of Incidence 
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Table A-9. Source Listing for Optical Analysis Program 




CQMiigN/M/ I PR . INTF, IBD.W(55.45) ,G_L. IZ.NCQAT, 1 BRANCH tPAPK 

00000010 


5 

COMM0N/REFL/TRANS(45) , C0Lt-(55j45) 

DIMENSION C0LLT(45) 

CONTINUE 

TOTWsQ ■ ■ — 

00000020 

00000030 

00000031 

00000032 ^ 



DO 7 1=1.45 
C0LLT( n=o 
DO 6 K=lj55 
COLLCK, I ) =0 

00000033 

00000034 

00000035 

00000036 

— 

6 

CONTINUE 

00000037 


7 ^ 

CONTINUE 

00000038 



CALL DATA IN 

00000050 



CALL MAIN 

00000060 



DO 20 1=1.45 

00000070 



DO 10 J=1,55 

00000071 



COLLTI I )=COLLT{ I )+COLLt J, I ) 

00000072 


10 

CONTINUE 

00000074 

— 


TOTW=TOTW+COLLT( I ) 

00000076 


20 

CONTINUE 

00000078 



DO 30 M=18j 90, 18 

00000080 



L=M-16 __ - - - 

00000081 



MM=M/2; LL=MM-8 

00000082 



WR1TE(6,90) (K-2,K,K=L,M,2) 

00000083 



WRI TEC 6,91 ) (COLLTCI ) , I=LL,MM) 

00000084 


30 

CONTINUE ‘ .. 

00000086 



WRlTEC6,9Z)T0Yvf 

00000088 


90 

FORMAT!/, 10C2X, 12, 12, IX)) 

00000090 


91 

FORMAT C9F8. 2) 

00000091 


92 

FORMAT! "TOTAL ENERGY COLLECTED= " , F9 . 2) 

00000092 

TT 



00000093 



PRINT "MAKING MORE RUNS? 1=YES, 0=N0“ 

00000110 



READ, I AGAIN 

000001 1 1 



IFCIAGAIN .EQ. 1 )G0 TO 5 

00000112 



STOF"' ■ 

00000113 



END 

00000114 

C 



00000130 ‘ 



SUBROUTINE MAIN 

00000140 




000001 50 

c 



00000160 



IBRANCH=1 

00000170 



INTF=1 

00000180 



IBD=1 

00000190 



IPR=1 

00000200 



IZ=0 

00000210 



CALL ONEPASS 

00000220 


10 

CONTINUE 

00000230 



IFdBRANCH .NE. 1 )G0 TO 40 

00000240 



IBD=1 

00000250 



INTF=2 

00000260 



■CALL TW6PA5S — — 

00000270 


20 

CONTINUE 

00000280 



IFCIBRANCH .NE. DGO TO 30 

00000290 



IBD=2 

00000300 



lNih=2 




CALL TWOPASS 

00000320 



IBD=3 

00000330 



CALL TWOPASS 

00000340 



iBRANUTTSU 




GO TO 20 

00000360 
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c 

c 

30 CgNTINUE : rr— — 

IPR=12 
12 = 9 
IBD=4 

Lb!IE=3 : ... — — 

CALL TWOPASS 
IZ=0 ^ 

IBD=5 

INTFsg : — 

CALL TW0PASS 
IBRANCH=0 
GO TO 10 

C : — 

40 CONTINUE 
IPR'=5 
IZ=25 

1 BRANCH=-1 : 

IBD = 6 
INTF=3 

CALL TWOPASS 

12 = 0. - — 

50 CONTINUE 

1F( I branch .NE. DGO to 60 
IBD = 7 

• 1NTF=2 ■ 

CALL TWOPASS 
1BD = 8 
INTF = 2^ 

. CALL TWOPASS » : _ ^ — — 

I BRANCH-0 
GO TO 50 
C 

60 CONTINUE ' r . — 

IPR=37 
12 = 9 
1BD = 9 

►INTF = 3 

‘ CALL TWOPASS 
IZ = 0 
1BD=10 

INTF = 2 - . . 

CALL'TWdPA"SS 

RETURN 

END 

.C ^ ^ - 

0 — ^ 

C 

SUBROUTINE ONEPASS 

CALL REFLCOL 

CALL REFRAcT"' — 

CALL SCATTER 
CALL ABSORB 

. . RETURN 

END 

C 

c 

C : - 

SUBROUTINE TWOPASS 

COMMON/M/IPRj INTF, I BD, W ( 55, 45 ) , GL, IZ,NC0AT, I BRANCH, PACK 


00000370 
00000380 
00000390 -.. 
00000400 
00000410 
00000420 
.0Qnp04.aQ_ 
00000440 
00000450 
00000460 
00000470 
00000460 
00000490 
00000500 
0000051.0 .. 
00000520 
00000530 
00000540 
00000550 
00000560 
00000570 
00000580 
00000590 
00000600 
00000610 
00000620 
00.000630. 
00000640 
00000650 
00000660 
0Q000S70 
00000675 
00000680 
00000690 
QQ000_7_0J1^ 
0000071 0 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000780 
00000790 
00000800 
00000610 
-00000820 
00000830“ 
00000840 
00000850 
00000860 
00000870 
00000880 
0000089Q 
00000900 
00000910 
00000920 
00000930 
00000940 
00000950 
00000960 
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c 


CALL DIVIDE 
CALL REFLCOL 

00000970 

00000980 

00000990 



GO TO (lOj 10j20j 10,20, lOj 10, 20j 10j20)j IBD 

00001 000 


10 

CONTINUE 

00001010 



CALL SCATTER 

00001020 



CALL ABSORB -- — — 

00001030 



INTF=4 

00001040 



CALL REFLCOL 

00001050 



CALL SCATTER 

00001060 



CALL ABSORB — 

00001070 


20 

CONTINUE 

00001080 



RETURN 

00001090 



END 

00001100 

c 



■ ■ 00_00JJ1J_Q 

c 



00001 120 

c 


. 

00001130 



SUBROUTINE REFLCOL 

00001 140 



COMMON/REFR/RF(S) .REFF 

00001150 



COMMON/REFL/TRANS ( 45 ) , COLL ( 55, 45) 

00001 160 



COMMON/M/ 1 PR, INTF, I BD, W(55, 45) ,GL, IZ, NCOAT, I BRANCH, PACK 

00001170 



DIMENSION E(45),G(45) 

00001 180 

c 



00001 190 



DTR=. 01 74533 

00001200 



GO TO (10,20,20, 10) j INTF 

00001210 ' 


10 

CONTINUE 

00001220 



IFCNCOAT .EQ. 0)GO TO 20 

00001230 



DO 11 1=1, 4S 

00001240 - 



G(I)=W( IPR, I ) 

00001250 



E(n-G( I )*TRANS(I ) 

00001260 


11 

CONTINUE 

00001270 




00001280 



DO 12 1=1,45 

00001290 



E(n=G( I)-E(I ) 

00001300 


12 

CONTINUE 

00001310 



GO Td 60 

00001320 

c 



00001330 


20 

CONTINUE 

00001340 



DO 22 1=1,45 

00001344 


22' 

G( I )=WUPR, I ) 

00001.345 



GO TO (24,25,26,27), INTF 

00001350 


24 

RFRN=RF(1 ) 

00001360 



RFRM=RF(2+NC0AT) 

00001370 



GO TO, 28 

00001380 


25 

RFRN=RF(2+NC0AT) 

00001390 



RFRM=RF(3+NC0AT) 

00001400 



GO TO 28 

00001410 


26 

RFRN=RF(2+NC0AT) ' ' 

00001420 



RFRM=RF(4+NC0AT) 

00001430 



GO TO 28 

00001440 


27 

RFRN=RF(2+NC0AT) 

* 00001450 




00001460 


28 

CONTINUE 

00001470 

c 



00001480 



DO 50 N=1,45 

00001490 



i=N+N-r ■ . 

COCCI 5oC 



R= 1 - C C RFRN/RFRM ) *S I N ( I *DTR ) ) * *2 

00001510 



IF(R .GT. 0)G0 TO 30 

00001520 



TR = 0 

00001530 



GO" rb”’40 



30 

CONTINUE 

00001550 
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m A _ 


ORtGiNAL PAGE jS 
OF POOR QUALITY 


c 


& 


RP= ( COS ( I *DTR ) - ( RFRN/RFRM ) *SQRT ( R) ) / ( COS ( I *OTR ) + 
(RFRN/RFRM)*SQRT(R)) 

BS=f nfisf t >!nTR1»PFBN/RF R M-snRT(R) ) / f COS ft ^PITRl *RFffN i ^RFRMt 


& SQRT(R)) 


TR=1 - (RP**2+RS**2)/2 


IFCTR .LE. 0)TR=0 

an CONTTNUE — 

IFCINTF .EQ. n SO TO 45 

GO TO (41 , 41 ,41 ,42j41 ,42,41 ,41, 42j41 ) j !BD 


41 TR=1-TR 

GO TO 45 

42 TR=1 -TR+REFF*TR**2 
45 E(N) =G(N)*TR 

50 CONTINUE 

60 CONTI NU.E : 


IPR=1PR+1 
DO 70 1=1,45 
MfIPR. n=E(l) 
70 CONTINUE 


C 

IFUNTF .EQ. DGO TO 85 

GO to ( ftfl 80 . 80 . 85. 80 - 85. 60, BO . 85. 801 ■ IBP 

80 CONTINUE 

DO 82 1=1,45 

ANGI = ( I +I -1 ) * . 0174533 

COLLf IPR. n = (Wf IPR-1 ■ n -M(lPR. m*COS<ANG-U 

82 CONTINUE 
85 RETURN 
END 

C ; : — : 

C 

c 

SUBROUTINE SCATTER 

COMMON/S/ SCAT(45,45,4) 

COMMON/M/IPR, INTF, IBD,W(55,45) ,GL, I2,NC0AT, 1BRANCH,PACK 
DIMENSION E(45),G(45) 

DO 10 1=1,45 

G ( I ) =W ( I PR, I ) • — 

E( I )=0 
10 CONTINUE 
C 

DO 30 1=1,45 i: — 

■ DO 20 J=1,45 

E(J)=E(J)+G(I )*SCAT( I , J, INTF) 

20 CONTINUE 

• 30 CONTINUE - 

1PR=1PR+1 
DO- 40 1=1,45 
W(IPR,1)=G( I )■••£( I) 

40 CONTINUE ■ — 

•RETURN 
. END 


C 

■ C - . - , . ... 

c 

. SUBROUTINE ABSORB 
COMMON/ABS/ABSF 

■ ■ COMMON/M/ 1 PR , I NT F. 1 BD ■ W(55 ■ 45) . 6L. IZ.NCOAT , 1 B RANCH. PACK 
DIMENSION E(45) 

DO 10 N=1,45 


00001560 

00001570 

■00001580 

00001590 

00001600 

00001605 

onnnifim 

00001820 

00001624 

00001630 

nnnni6GP_ 

00001634 

00001640 

00001650 

nnnmsBO- 

00001670 
00001680 
00001690 
■jaQnoi7QO- 
00001710 
00001720 
00001725 
00001730. 
00001740 
00001750 
00001755 
00001760 
00001770 
00001780 
00001790 
00001 800_ 
00001 810 
00001820 
00001830 
00001840 
00001850 
00001860 
00001 870 
00001880 
00001885 
00001 890 
00001900 
00001910 
00001 920 
00001930 
00001940 
00001950 
00001960 
00001970 
00001980 
00001990 
00002000 
00002010 
00002020 
00002030 
00002040 
00002050 
00002060 
00002070 
00002080 
00002090 
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E(N)nWUPR.N) 

l=N+N-V 

3sQL/C0SC I 0174503^ - ^ — 

00002100 

00002110 

000021^ 


10 

E(N)=E)<P(-S*ABSF)*E(N) 

CONTINUE 

IPR=IPR+1 

DO 20 1 = 1.45 — . 

00002130 
00002140 
000021 SO 
00002160 

n 

20 

WCIPH^ i )=EU ) 

RETURN 

END 

00002170 

00002180 

000021.90 

00002200 

c 



00002210 

n 



00002220 



SUBROUTINE REFRACT 

00002230 



COT1MON/.REFR/RF<5),REFF .. - 

00002240^ 



COMMON/M/lPRx INTF, IBDj W( 55, 45) , QL, IZ,NCOAT, IBRANCH.PACK 

OCi002250 



DIMENSION G(45),E'(45) 

00002260 

Q 



00002270 



IFMNTF .NE, 1 ). 80 TO 40 

00002260 




J=1 

00002290 



DO 5 1=1,45 

00002292 



E( 1 ) =0 

00002294 


5 ^ 

8(1 )=W( IPR, 1 ) ' - 

00002296 



DTR='.'i5T'74533 

00002300 



RFRN=RFn ) 

00002310 



RFRM«RF(2+NCDAT) 

00002320 



DO 20 N=1., 45 

00002330 



I=N+N-1 ' 

00002350- 



A=ARSIN(SIN( ( 1-1) *DTR)*RFRN/RFRM)/DTR 

00002360 ' 



C=ARSIN(SIN(( 1+1 >*DTR)*RFRN/RFRM)/DTR 

00002370 


11 

CONTINUE . - 

00002380 

— 


iF(A -LT. J+1 . )GO TO 12 

00002390 



J = J + 2 

00002400 



GO TO n 

00002410 


12 

CONTINUE. 

' 00002420 



IF(C .GT. J-1.)GO TO 13 

00002430 



J=J’“2 

0000244C- 



GO TO It 

00002450 


10 

CONTINUE . - 

00002460 



K=( J+1 )/2 

00002470 



IF(C ,LT. J+1 . )G0 TO 14 

00002480 



E(K+1 >=E(K+1 )+(C-J-1 )s<G(N)/(C-A) 

00002490 



E(K)=E(K)+(J+1“A)*G(N)/(C-A) 

00002500 



GO TO 20 

00002510 


14 

contiiQue 

0000252Ci 



IF(A .GT. J-1 . )G0 TO 15 

00002530 



!E(K)=E(K)+(C-J+1 )*G(N)/(C-A) 

00002540 




00002550 



GO TO 20 

00002560 


IS 

E(K)=E(K)+G(N) 

00002570 


20 

CONTINUE 

00002580 

C“ 






IPR=IPR+1 

00002600 



DO 30 K=l,45 

0000281 0 



W(1PR,K)=E(K) 

00002620 

■ 

-“-0ior 

“UONUNUb 



40 

RETURN 

00002640 



END 

00002650 

c 



00002660 

xr 




0 



00002680 
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SUBROUTINE DIVIDE 
COmON/REFR/RF ( 5) , REFF 

aeiMWQN/MZlFJR . 1 NTF , T BD , W f 55^451, f5L^JL7^,ND0AT^RRANDIi,_FAnK ‘ . 

00002690 
00002705 
nnon5>7.i o 


DIMENSION £(45)^6(45), H(4S) 

00002720 

c 


00002730 


m 10 1=1,45 

00002740 


GLai^wcrpR.i) 



1ANG=I 

00002755 


Em=em*CELLBDF(IBD,lANG,GL, NCOAT , PACK , ,RF ) 

00002760 


H<n=em-E(i) 

00002770 


rOCONilNUE . 

1 . a— 


IPR=!PR+1+IZ 

00002790 


IF( [BRANCH .EQ. DGO TO 30 

00002800 


DO 20 1=1,45 

00002610 


wupR.nsHm > . 



20 CaNTINUE 

00002830 


eo ra so 

00002840 


30 40 'I'=l ,45 

00002650 


w(ipR.n=:Ecn . . 



40 CONTINUE 

’ 00002870 


50 RETURN' 

00002880 


END ' 

00002890 

c 


■■ MM a—— 

c 


00002910 

c 


00002920 

* 

SUBROUTINE DATA IN 

00002930 


COMMON/REFL/TRANS ( 45 ) , COLL( 85 . 45 ) 

- 00002940 


COMMON/ABS/ABSF 

00002950 


COMMON/REFR/RF ( 5 ) , REFF 

00002960 


COMM0N/M/IPR, INTF,4BDj W(55,45X,GL, IZ, NCOAT. IBRANCH. PACK 

00002980 


COMMON/S/SCAT C 45 , 45 , 4 ) 

00002990 


■ DIMENSION FSCAT(S), SOLAR (45) 

00003000 


DATA FSCAT/ . 155, .17,. 035, . 045, , 005/ 

00003010 



00003020 

0 

INSOLATION FOR PHOENIX 1953 

00003030 


DATA sOLar716.7, 20 . 1 , 1 9. 3, 1 9 . 3, 20. 5, 26 . 2, 43 . 1 , 1 30. 2, 90 . 9, 

00003040 


a 101. 1,54. 1,29. 4, 27. 4, 30. ,144, 8,173. ,101, 3,68. ,35. ,.21. 6, - 

00003050 


--a 24. 2,. 106 ..2, 176..,7„24.6,21 .4. 17. 9, 20. t , 26. 5, 38. 1 , 103. 9, 1 9. 3, 

00003060 


a 1 1‘. 4,9. 1 , 1 1 .7, 4.5, 12.5,33.3, 12. 1 .2.6, . 8, .4, . 0, 1 . 7, 1 .7, . 5/ 

' 00003070 


- 

00003080 

c 

SCATTERING FUNCTION LOADING 

00003090 

c 


00003100 


‘DO 10 1=1,45 

000031 TO 


SCATC I, I, 1 )=0. 

000031 20 


SCATC I , r,2)=-*94 

00003130 


SCATCI , I, 3) =-..94 

. 00003140 

_ 

'SCATd, i,4) = -,oa 

00003150 


10 CONTINUE 

00003160 


DO 20 1=1,44 

000031 70 


. SCAin+l, t,4) = ,04 .. 

00003180 


SCATCd I+l .4)=. 04 

000031.90 


20 CONTINUE ' 

00003200 


DO 40 J=l,5 

0.0003210' 


K=45-J 

00003220 


FK=FSCATCJ) 



• DO 30 1=1, K . . ' . 

00003240 


'FK1=FK, 

00003244 


•IF(1 '.LE. 5. .AND. I +U .GE. 2*1 )FK1 =2*FK 

. '00003245 


SCATCI+J, I ,25=FK ’ ' 

00003250 


SCAT< 1 , l+J,2)=FK1 ' '' 

00003260 

- ^ 

SCAT< I+J, i ;31s:FK ’ 

, 00003270 
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c 

c 


c 


c 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

0 

c 


SCATCI, I+J,3)=FK1 
30 CONTINUE 

40 CONTINUE 

INSOLATION DATA LOADING 
DO 45 J=1,45 


00003280 

00003290 

QJO.OQ33QO 

00003310 

00003320 

00003330 


_ANGIg(J-i-J-1 )«.Q174S3 


-.OLoaogaas.. 


45 • W( 1, J)=SOLAR( J)/COSCANGI) 

PRINT "NO. OF COATING LAYER, 1 OR 0" 


00003340 

00003350 

00003360 

■00003370 


IFCNCOAT .EQ. 0)G0 TO 60 

PRINT "COATING INDEX OF REFRACTION" 
_READ.RF(2) 


00003380 

00003390 

00003400 

-QQD034,I0 


PRINT "COATING THICKNESS & WAVELENGTH (MUST HAVE SAME UNITS)" 
READ,COATL, WAVEL 
60 CONTINUE 
SF Q ig. l , 


00003420 

00003450 

00003460 

-QQQQ347.Q, 


PRINT "INDEXES OF REFRACTION OF GLASS COVER, BONDING ADHESIVE 
&ENCALSULANT" 

READ, (RFC J+NCOAT) , J=2,4) 

PRINT "MODULE PACKING FACTOR (BETWEEN ,512 8 .614)" 


e 0003480 
00003490 
00003510 
00003520 


READ, PACK 00003530 
PRINT "GLASS THI CKNESSO. 2, 4 OR 4.8MM), BULK ABSORP. COEFF . (1 /MM) " 3540 
READ, GL, ABSF 00003560 
PRINT "REFLECTANCE OF ENCALSULANT" 00003562 


?^EAD, REFF 

CALCULATION OF TRANSMISSIVITY FOR ONE LAYER COATING 


00003564 

00003570 

00003580 

00003590 


IFCNCGAT .EQ. 0)GO TO 80 
DTR=. 017453 
DO 70 N«l,45 

I=N-^N-1 


00003600 

00003610 

00003620 

00003630 


INCIDENCE ANGLES IN MEDIA FROM SNELL'S LAW 


00003640 

00003650 

00003660 


ANG1-UDTR 

ANG2=ARSIN( (RFC 1 )/RF( 2) ) J«cSI NCANGl ) ) 
ANG3=ARSI N( (RFC 1 )/RF(3) )*SIN(ANG1 ) ) 


EFFECTIVE INDEXES FOR NON -NORMAL INCIDENCE 

RFISgRFCl )*C0SCANG1 ) 

RF2S=RF ( 2 ) COS ( ANG2 ) ^ 

RF3S = RF ( 3 ) * COS ( ANG3 ) 

RF1P=RF(1 )/C0SCANGl ) 

RF2P = RF C 2 ) /COS C ANG2 ) 

RF3P=:RF(3) /COSC ANG3) 

PHASE ANGLE BETWEEN AMPLITUDE VECTORS 

PHS=(4^0. 1416*CQATL/WAVEL)»RF2S 

AMPLITUDE REFLECTANCE S & P COMPONENTS 
R1 Sg ( RF2S -RFl S ) / ( RF2S+RF1 S ) 

R2S= C RF3S-RF2S ) / ( RF3S+RF2S ) 

R TP= ( RF2P - RFl P ) / ( RF2P+RF1 P ) 

R2P= ( RF3P - RF2P ) / ( RF3P+RF2P ) 


J>S>.QQ^^.7Q 

00003680 

00003690 

00003700 

00003710 

00003720 

00003730 

00003740 

00003750 

00003760 

00003770 

00003780 

00003790 

00003800 

00003810 

00003820 

00003830 

00003840 

00003850 

00003860 

00003870 

00003880 

00003890 
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c 

C INTENSITY REFLECTIVITIES 



RP = R1PJK){^2+R2P«*2+2*R1P5KR2P*C0S(PHSD 
C 

C AVERAGE FOR P0L0RI2ED LIGHT 

R=(RS+RP)/2. - - - 

TRANS(N)=1-R 
70 CONTINUE 
80 CONTINUE 

RETURN - 

END 

FUNCTION CELLBDFC IBD, I , GL^ NCOAT, PACK^ RF) 

DIMENSION RSHIFTanO) , RSHIFT1 ( 1 0) , FPACK(4 ) , RFC 5) 

DIMENSION DD(10,4) .DND(10.4) >NPP(10s4) .NPND(10,4) .DDDC10.4) , 
SDNDDC 1 0 , 4 ), NDDDC lOj 4 ) , NDNDDC 10^4) 

REAL NDDjNDND,NDDD,NDNDD 
C 

DATA RSHIFT1/0> , > 01 , . 02. . 03, . 04 , « 05, , 06, , 07, . 08, . 1 / 

DATA RSHIFT2/. 0, .02, . 04 , . 06, . 08, , 1 ^ . 1 2, . 1 6, . 1 8, . 2/ 

DATA FPACK/.512, .658, .782, .814/ 

C 

DATA DD/1 . . .994. .98, .965, .951 . .936. .922. ,907. .893, .868. 

81 . . . 994. . 98. . 965. . 951 . . 936. . 922. . 907. . 893. . 868, 

81 . . . 994. . 98. . 965. . 951 . . 936. . 922. . 907. . 893. . 868, 

81 . „ . 994, . 98, . 965, . 951 , . 936, . 922, . 907, . 693, . 868/ 

*C 

DATA DDD/1 . , . 992, . 965, . 934, . 895, . 846, . 8, . 705, . 66, . 63, 

81 . . . 992. . 965. . 934. . 895. . 646. . 8. . 705. . 66. . 63, 

8T. , . 992, . 965, . 934, . 9, . 873,, . 845, . 79, . 763, . 736, 

81 . , .992, .965. .934, .9, .873, .845, .79, .763, .736/ 

C 

DATA DND/0., .006, .021, .042, .061, .082, .1, .117, .134, .15, 

80.. .008. .027. .048. .071. .094. .116. .139. .161. .183, 

8.0. .008. . 027. .048. .065. . 081 . . 1 . . 1 1 7. . 1 34 . . 1 5 , 

8.0. , 006, . 02, . 036,'. 051 , . 064, . 075, . 09, . 1 02, . 1.1 5/ 

C 

DATA DNDD/0. ,0.,0,,0.,0.,0.,0.,0.,0.,0., 

80. . 0. . 0. . 0. . . 001 . .0054. . 0156. .0156. . 0156. . 0156, 

so; ,0. , .0012,7006, .0122, .02, .0296, .0296, .0296, .0296, 

80. . . 0014. . 0048. . 0096. .0158. .023. . 032. . 032. . 032. . 032, 

C 

DATA NDD/0. , ,008, .025, .046, .07, .096, . 124, ■ 124, . 124, ■ 124, 

SO. ,.01 6-, .044, .08, .12, .164, .200, .208, .208, .208, 

80. . .025. .076. . 133. . 19. .24. .29. .29. .29. .29, 

80 . . . 025. . 076. . 1 33. . 1 9. . 24 , . 29, . 29, . 29, . 29/ 

C 

DATA NbOb/O . , . 01 2, . 03377^^^ * 081 , . 105, . 105, . 105, . 1 05, . 105, 

80. . . 012. .042. . 079. . 1 1 6. . 1 53. . 1 53. . 1 53. . 153. . 153, 

8. 0. . 012. .06. . 1 18. . 172. .24. .24. .24. .24. .24, 

80. . . 012. . 06. . 1 16. . 172. .24. .24. .24. .24. .24, 

DATA NDND/1 . , . 99, .967, .94, .912, .88, .855, .855, .855, .855, 

81 . . .986. .945. .892. .84. .798. .767. .767. .767. .767, 

81 . . . 962. . 88. . 79. . 71 . . 64. . 585. . 585. . 585. . 585, 

. 97; : 837, . 81 , .725, . 656: 

C 

DATA NDNDD/0. , .01 , . 020, . 051 , . 076, . 1 01 , . 1 01 , . 1 01 , . 1 01 , . 1 01 , 

80. . . 018. . 051 . .078. . 098. . 108. . 108. .108. .108. .108, 

80. , .018, . ‘ Obv; ; ‘ (r8,V054>;Vr D SV7TO5,7Tggr "^ ^ 

80 . . . 025. . 058. . 082. . 095. . 1 05. . 1 05. . 1 05. . 1 05. . 1 05, 


00003900 

00003910 

0Q0Q392Q. 

00003930 

00003940 

00003950 

00003970 

00003960 

00003990 

00004200 

00004210 

00004250 

00004260 

00004270 

00004280 

00004283 

00004285 

00004290 

00004300 

00004310 

00004315 

00004320 

00004330 

00004340 

00004350 

00004355 

00004360 

00004370 

00004380 

00004390 

00004395 

00004400 

00004410 

00004420 

00004430 

00004435 

00004440 

00004450 

00004460 

00004470 

00004475 

00004480 

00004490 

00004500 

00004502 

00004505 

00004510 

00004520 

00004530 

00004540 

00004545 

00004550 

00004560 

00004570 

00004580 

00004585 

00004590 

00004600 

06004610 

00004620 
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Table A-9. Source Listing for Optical Analysis Program (Cont’d) 


= 017453 

RFMgRF/g) _ 

IF<NCOAT .EQ. 1)RFM=RF<3) 

ANGS=ARSIN(SIN(A1)/RFM) 

S=2*6LJKTAN(ANGS)/5.954S/aS.4 

gg TO t10.12,14.1S.18 .?O-gg.g4.gS.g8>.IBO 

10 CELLBDF=PACK 

eg TO 35 

12 CELUBDF=TNT2{3, PACK^ 10,4, RSHIFT1 , FPACK, DO, II , 12, 1 O) 

60 TO 35 

14 S=S*2 

CELLBDF=TNT2(2,PACK, 10,4, RSHIFT2, FPACK, ODD, 1 3, 14, 10) 
GO. TO 3S 

16 CELLBDF=TNT2(S.PACK. 10.4.RSHIFT1 . FPACK. DND, 15. 16. 10) 
GO TO 35 
18 S=S*2 

CELLBDF=TNT2(S,PACK, 10,4,RSH1FT2,FPACK,DNDD, 17, 18, 10) 

GO TO 35 

20 CELLBDFsl -PACK 
GO TO 35 

• 22 CELLBDF=TNT2(S,PACK, 10,4,RSHIFT1,FPACK,NDD, IS. no, 10) 
GO TO 35 

24 S=?S*2 - , ..w 

CELLBDF:5TNT2 (S, PACK, 10,4, RSHI FT2, FPACK, NDDD, JT , 10) 

80 TO 35 

26 CECLSDF«TNTaCS.PACK.10«4.RSHIFT1 .FPACK.NDND.J3. J4, 10) 
GO TO 35 
28 S^S*2 

CELL8DF=TNT2CS, PACK, 10,4, RSHIFT2, FPACK, NDNDD, J5, J6, 1 0) 

35 CONTINUE 

RETURN 

END 


00004625 

00004630 

..,00X10.4, 64 0., 

00004650 

00004660 

00004670 

.QQQD 468Q , 

00004690 

00004700 

00004710 

J20-Q0,4Z20.. 

00004730 

00004740 

00004750 

00004770 

00004780 

00004790 

00004800. 

00004610 

00004620 

00004830 

00004840 

00004850 

00004860 

00004870 

smnm^3.Q^ 

00004890 

00004900 

00004910 

QQ004920 

00004930 

00004940 
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